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ABSTRACT 

Wc present deep 1.2 mm continuum mapping of a 566 arcmin 2 area within the Lockman Hole 
North field, previously a target of the Spitzer Wide-area Infrared Extragalactic (SWIRE) survey 
and extremely deep 20 cm mapping with the Very Large Array, which we have obtained using 
the Max-Planck millimeter bolometer (MAMBO) array on the IRAM 30 m telescope. After 
filtering, our full map has an RMS sensitivity ranging from 0.45 to 1.5mJybeam -1 , with an 
average of 0.75 mJy beam -1 . Using the pixel flux distribution in a map made from our best data, 
we determine the shape, normalization, and approximate flux density cutoff for 1.2 mm number 
counts well below our nominal sensitivity and confusion limits. After validating our full dataset 
through comparison with this map, we successfully detect 41 1.2 mm sources with S/N > 4.0 and 
Si. 2 mm — 2 — 5 mJy. We use the most significant of these detections to directly determine the 
integral number counts down to 1.8 mJy, which are consistent with the results of the pixel flux 
distribution analysis. 93% of our 41 individual detections have 20 cm counterparts, 49% have 
Spitzer/MIPS 24 (im counterparts, and one may have a significant Chandra X-ray counterpart. 
We resolve ~ 3% of the cosmic infrared background (CIB) at 1.2 mm into significant detections, 
and directly estimate a 0.05 mJy faint-end cutoff for the counts that is consistent with the full 
intensity of the 1.2 mm CIB. The median redshift of our 17 detections with spectroscopic or 
robust photometric redshifts is z m odian = 2.3, and rises to ^median = 2.9 when we include rcdshifts 
estimated from the radio/far-infrared spectral index. By using a nearest neighbor and angular 
correlation function analysis, we find evidence that our S/N > 4.0 detections are clustered at the 
95% confidence level. 
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1. Introduction 

Over a decade ago, measur ement of the cos- 
mic infrared backg round (CIB: iPuget et al.lll996t 
Fixsen et al.l 119981 ) revealed that approximately 
half of all of the light in the universe emerges 
at far-infrared wavelengths due to reprocessing 
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by dust (e.g., iDole et al.l 120061 ). With the ad- 
vent of the Submilli meter Common- User Bolome- 
ter Array (SCUBA; iHolland et aD Il999h . it be- 
came clear around the same time that not all 
of this reprocessed emission originated in galaxy 
populations that could be easily detected with 
optical telescopes. Instead, SCUBA surveys at 
850 yum revealed the existence of a p opulation of 



bright submillimeter galaxies (SMGs: ISmail et ah 
1997t iBarger et alJll998t iHughes et al.lll998l) with 



faint or undetectable optical counterparts. Opti- 
cal spectroscop y of the least obscured m embers of 
the population ( Ivison et al. 1998 . 2000l). followed 



up by dete ctions of CO emission (jFraver et al 



19981 Il999t) . confirmed that SMGs were indeed 
a high-z population. Their faint X -ray coun- 
terparts ( Alexander et al. 20031 20051) as well as 
mid- infrared spectroscopy (IValiante et al. 20071 



Menendez-Delmestre et~al1l2007ll2009tlPope et"aT 
20081 ) indicated that they were not predom- 
inantly powered by accretion, but rather by 
star formation. These observations suggest that 
SMGs may play an important role in the cos- 
mic star formation history. However, the details 
of their star formation remain uncertain. Dif- 
ferent schools of thought exist about whether 
SMGs are major mergers ( e.g. , IConselice et al 



20031 iNaravanan et al.ll2009l l2010t ) or host bursts 
triggered by the inflow of intergalactic gas along 
filaments ( Dave et al.ll201oh . Likewise, galaxy evo- 
lution models that con sider SMGs have d isagreed 
on whether they are (IBaugh et al. 120051) or arc 



not ([Havward et al.l 120111) forming stars with an 
unusually top-heavy initial mass function (IMF). 

Understanding how SMGs fit into the over- 
all history of cosmic star formation has been 
impeded by their high obscuration in the op- 
tical and the coarse angular resolution of the 
(sub)millimctcr bolometer arrays used to detect 
them. An important advance came with the use 
of deep, high-resolution radio continuum maps 
with the Very Large Array (VLA) to localize 
SMGs on the basis of the far-IR/radio correla- 
tion. K eck followup of radio-preselected SMG 
samp les (jlvison et al.ll2002l : IChapman et al.ll2003 . 
20051 ) broke the logjam for determining SMGs' 
redshifts, allowing a dramatic expansion of SMG 
samples with CO-confirmed spectros copic red- 
shifts and spatially resolved maps (jNeri et al 
20031 : iGreve et al.l 120051 : iTacconi et al.ll2006l ) Al- 



though these developments have been important 
and exciting, there are a number of caveats on 
the current state of our knowledge. First, not all 
SMGs have counterparts in VLA maps of typ- 
ical depths, and not all that have counterparts 
yield optical spectroscopic redshifts; this has led 
to uncertainties in the population's overall red- 
shift distribution, especially at the high-z end. 
Second, we have come to realize that the 850 /im 
waveband at which the influential early SCUBA 
work was done does not give us a complete picture 
of all dusty galaxy populations at high rcdshift: 
observations at shorter/longer wavelengths pref- 
erentially detect populations with higher/lower 
dust temperatures and/or l ower/higher redshifts 
e.g.. IChapman et al.l 12004 IValiante et al.l 120071 



Magnelli et al.lbOloHChapman et al.l2010HMagdis et al" 
20101 ). Finally, with the exception of highly lensed 
systems, our direct knowledge is restricted to 
bright individual detections for which limited sen- 
sitivity and angular resolution (vs. confusion) 
are not fatal obstacles. Understanding the lower- 
Lir galaxy populations that produce the bulk of 
the CIB is challenging, and efforts to determine 
their properties are not al ways consistent wi th the 
CIB's normalization (e.g., Scott et al. 2010() . 

To fill in some of the gaps in our knowledge, 
we need to obtain deep mapping at multiple 
(sub)millimcter wavelengths, at high angular reso- 
lution, over a large area that has good multiwave- 
length coverage (and especially, very deep radio 
continuum data). This combination of properties 
can in principle allow us to (a) optimize the identi- 
fication of counterparts at other wavelengths, and 
therefore the determination of redshifts and the 
assessment of SMGs' detailed evolutionary states; 
(b) defeat cosmic variance and optimize sensitiv- 
ity to clustering and large-scale structure, a key 
first step in connecting the properties of SMGs to 
the properties of the dark matter halos that host 
them; and (c) strongly constrain the parameters of 
SMG number counts down to faint flux densities, 
so that we can accurately compare the census of 
obscured star formation to the constraint of the 
CIB. 

In this paper, we report 1.2 mm observations at 
11" resolution of a unique deep field that sat- 
isfies many of these needs. Our map is both 
larger and more sensitive t han previous deep 
maps made at 1.2 mm (e.g., IGreve et al.l l2004t 
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Bertoldi et al.ll2007t iGreve et al.ll2008l ). and com- 
pared to recent work done with other instruments, 
displays a competitive combination of attributes. 
Surveys made at 1.1mm using ASTE/AzTEC 
have mapped equal or wider fields to a greater 
depth , but with much lower spatial resolutio n 



(e.g., IScott et all 120101 lHatsukade et al.l l201lh 



Maps using JCMT/AzTEC and APEX/LABOCA 
generally achieve wider fields or greater depths, 
but not bot h, and also have coarser spatial reso- 
lution (e.g. , IScott et al.l 120081: IPerera et al.l 12008 ; 



Weifi et all 120091 : lAustermann et al.l l2010h . The 
powerful combination of resolution, depth, and ex- 
tent of our MAMBO map, together with the rich 
supplementary data available for our target field, 
make it a powerful tool for studying the properties 
of SMGs. 

Our map lies within the Lockman Hole North 
(LHN) field, centered at 10:46:00 and +59:01:00 
(J2000), which was one of the targets of the Spitzer 
Wid e- Area Infrared Extr agalactic (SWIRE) sur- 
vey (|Lonsdale et al.ll2003l) . In addition to SWIRE 
coverage in all four IRAC and all three MIPS 
bands, the LHN has been the subject of extremely 
deep 20 cm continuum mapping with the Ver y 
Large Array (VLA) by I Owen fc Morrison! (|2008f ). 
who produced an ABCD configuration synthe- 
sis image with a central ler RMS sensitivity of 
2.7 /uJy. These data have bee n further supple- 
mente d by 90 cm VLA mapping (lOwen fc Morrison 
2009a!) . 50 cm GMR T mapping (jFiolet et al.ll2009t 



Owen et al 



(Owen et al 



201 lal). deep Spitzer 24 /im imaging 
2011b), optical sp ectroscopy with 



WIYN (|Owen fc Morrison! l2009al ). X-ray imaging 
from the Chan dra /SWIR E survey (|Polletta et al 



20061: IWilkes et al.ll2009l) . and determination of 



photometric rcdshi fts from multicolor opt ical and 



ncar-IR imaging (Str azzullo et al.l I2010I ). The 
LHN is also one of the targets of the Her- 
schel Multi-tiered Extragalactic Survey (HerMES; 
Oliver et al1l2010h . 



Observations 

We used the 117-element Max-Planck millime- 



ter bolometer (MAMBO) array (jKrevsa et al 



19981 ) at the IRAM 30 m telescope to obtain on- 
thc-fly mapping of the LHN at an effective wave- 
length of 1.2 mm. Our observations were obtained 
over the course of five semesters from 2006 through 



2010 (Table [l}. Due to telescope control software 
problems during the first two semesters - an error 
in computing corrections for atmospheric refrac- 
tion, which undermined the quality of the point- 
ing during the first semester, and tracking jitter 
that undermined map reconstruction during both 
semesters - we have restricted our initial analy- 
sis (e.g., P to a "best" map that includes only 
the data from our final three semesters of obser- 
vations. We have also constructed a "full" map 
using data from all five semesters, whose reliabil- 
ity we can validate based on comparison with the 
"best" map ( £I5.2[) . and which we therefore use for 
our analysis of the bright source population in the 
field. 101.3 hr and 192.5 hr of data were combined 
to produce the "best" and "full" maps with ar- 
eas of 514arcmin 2 and 566arcmin 2 , and average 
depths of 0.90 mJy beam -1 and 0.75 mJy beam -1 , 
respectively. 

During all five semesters, our MAMBO data 
were acquired during the weeks that IRAM dedi- 
cates to pooled observations of multiple bolometer 
projects. As a result, most of our data were taken 
with 250 GHz zenith opacity < 0.3, low sky noise, 
and essentially no cloud cover. Observations were 
limited to LST ranges when the LHN had eleva- 
tion > 40° (to minimize opacity corrections and 
pointing anomalies) and < 80° (to avoid slewing 
errors and sudden accelerations of the MAMBO 
array) . We built up coverage of our field by mak- 
ing many small ~ 300" x 320" (azimuth x ele- 
vation) maps, each of which required 41 minutes 
to complete. The pointing centers of these small 
maps were arranged in a rectangular grid with 2' 
spacings between map centers. 

We planned the observations so that before 
(and usually after) each 41-minute map, the tele- 
scope was pointed on a nearby quasar- typically 
J1033+609 at a distance of 1.9° from field center, 
with Si. 2mm ~ 0.3 Jy. Every 1-2 maps, a sky- 
dip was performed to measure the zenith opacity, 
confirm that weather conditions remained good, 
and allow the interpolation of opacity corrections 
during the maps themselves. Standard flux cal- 
ibrators were observed roughly every four hours; 
these coincided with large slews, to force resets 
of the telescope's inclinometers. Focus measure- 
ments on bright quasars or planets were made at 
the beginning of each observing session as well as 
after sunrise and sunset. We required all of these 
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calibrations so as to minimize pointing errors and 
anomalous refraction, as is important for the de- 
tection of faint point sources in a wide-area map. 

The IRAM 30 m uses a chopping secondary mir- 
ror to subtract low-frequency sky noise from on- 
the-fly MAMBO maps. This chopping gives the 
telescope an effective double-beam point spread 
function (PSF) on the sky, with one positive and 
one negative beam separated in azimuth and sym- 
metric about the nominal pointing center. Dur- 
ing shift-and-add (SAA) reconstruction (see Sj3|, 
negative-beam data are inverted and combined 
(for a given sky position) with their positive-beam 
counterparts. The end result is a triple-beam pat- 
tern that is a well-defined function of position 
for any single observation: two negative sidelobcs 
bracket a positive beam in azimuth. The SAA 
algorithm thus conserves the mean flux of the ob- 
servations, in that the negative sidelobes together 
contain as much integrated flux as the positive 
beam. When many observations are combined into 
a single mosaicked image, the effective PSF is a su- 
perposition of many triple-beam patterns that can 
vary substantially with position. 

Since a given sky position within the LHN usu- 
ally falls within several of our small maps, varying 
the chop throw and scan direction tends to scatter 
the negative flux into a uniform annulus around 
the positive central Gaussian, reducing its peak 
intensity and minimizing its deleterious effects on 
the fluxes of nearby pixels. We therefore (a) used 
different chop throws for alternating columns in 
our grid of pointing centers, and (b) recorded the 
scan direction of each map in equatorial coordi- 
nates, so that observers (to the best of their abil- 
ity) could observe new grid positions at LSTs such 
that scan directions would not match those of (ex- 
isting) adjacent maps. In the central regions of our 
final maps (i.e., where we extract sources), peak 
intensities of the negative sidelobes reach only 5% 
of the peak positive flux thanks to our adoption of 
these strategies. 

During our first two semesters of observations, 
we obtained maps of 88 grid positions out of 100 
in a 10 x 10 grid centered on the LHN field cen- 
ter stated in gTJ During our last three semesters, 
which contribute to the "best" map, we observed 
97 positions of 99 in a 9 x 11 grid, extending 2' 
farther east but 2' less far south than the origi- 
nal grid, as well as two additional grid locations 



in the southeast corner. Between the first and sec- 
ond semesters, we swapped which sets of pointing 
centers were observed with which chop throws (36" 
and 48" throws vs. 42" and 36" throws for alter- 
nating columns). Due to the differences in spatial 
coverage and weather conditions during the obser- 
vations, the areas where the "best" map and the 
"full" map are respectively deepest overlap but do 
not match perfectly (see Figure [TJ. 

3. Data reduction 
3.1. Signal maps 

The raw bolometer time stream data were re- 
duced using Robert Zylka 's MOPSIC0 pipeline, 
which is distributed in parallel with IRAM's 
GILDAS package. MOPSIC is the standard 
package for r educing deep MAMBO on-the-fly 
maps (see e.g., Greye et "all 20041: Voss et al. 2006t 



iBertoldi et al]|2007t IGreve et al.ll2008l) . We now 
briefly outline the steps of the MO PSIC reduc 
tion pipeline; for further details see IGreve et al 



(|2004| ). The pipeline removes spikes in the time 
stream data stronger than 5x the instantaneous 
bolometer RMS noise. It also subtracts a third- 
order polynomial baseline in time and performs 
correlated signal filtering on the bolometer time 
streams to identify and remove foreground at- 
mospheric emission that affects many bolometers 
simultaneously. Each bolometer is correlated with 
an annulus of neighboring bolometers within a 
150" radius, and the average signal of the twelve 
most highly correlated bolometers is subtracted 
away. The filtered time streams are then binned 
into 3.5" x 3.5" pixels, and a signal map is re- 
constructed using the SAA algorithm. The indi- 
vidual signal maps are combined into a mosaic 
image by averaging the map flux density at each 
pixel weighted by the local inverse variance. Our 
"optimally filtered" signal map (Figure [2]) was cre- 
ated by applying a PSF-matched filter to the final 
mosaic image f £|5.1[) . 

In addition to the signal map, the MOPSIC 
pipeline also produces a weight map that is locally 
proportional to the inverse variance in the signal 
map. By enforcing that the Gaussian distribution 
of the S/N map pixel distribution has a standard 
deviation of unity, we normalize the weight map 



1 scc http: //www. iram. es/IRAMES/mainWiki/CookbookMopsic 
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so that it can be used to find the local RMS noise, 
a = 1/VW, in the image (see Figure [1]). Using 
the weight map as a guide to find the local RMS 
noise for a detection is more robust than using the 
nearby pixels themselves, because locally the pix- 
els are affected by the negative residual sidelobes 
of SAA reconstruction, as well as those of other 
bright nearby sources. 

3.2. Noise maps 

Because of the telescope's effective triple-beam 
PSF, each source in the field injects negative as 
well as positive flux into the map. To generate re- 
alizations of source- free maps, hereafter referred to 
as "noise maps," we removed the negative and pos- 
itive flux from undetected as well as bright sources 
using two different techniques. We go on to use the 
different results for different purposes. 

We constructed the first type of noise map 
with a techniqu e common in MAMBO data anal 
ysis (see, e.g., 



Greve et all [2004; Bcrtold i et al 



2007tlGreve et al.ll2008l ). using the data reduction 
pipeline to scramble the known locations of the 
bolometers within the image plane. During recon- 
struction of the time stream data, this has the ef- 
fect of smearing the flux from any one source into 
an area on the sky of approximately 200arcmin 2 , 
reducing the intensity of the source's peak flux 
contributions by a factor of ~ 10 3 and making the 
peak flux contribution from our strongest sources 
~ 200 times fainter than the RMS noise. Because 
the telescope's chopping ensures that the mean of 
the map is zero, there is no residual baseline in- 
crease as the negative flux contributions arc iden- 
tically smoothed. These "shuffled noise maps" arc 
simple to construct, but it is cumbersome to pro- 
duce large numbers of them since each requires 
a full reduction of the data using the MOPSIC 
pipeline. Therefore, we use the shuffled noise maps 
only to estimate the noise of our "full" data dur- 
ing source extraction f £|5.ip as well as in the Monte 
Carlo simulation of completeness ( 35.5[l . 

We needed to develop a different technique for 
creating noise maps in order to quickly gener- 
ate thousands of independent noise realizations 
of chopped data for our P(D) analysis. For this 
we subtracted subsets of the data that are "jack- 
knifed" in the sense that we remove fractions of 
the original data first. One full image of our field 
is created using the data from only one bolome- 



ter in the array at a time. All bolometers other 
than the one of interest are masked away after the 
correlated signal filter is applied, so the data still 
receive the benefit of correlated sky noise subtrac- 
tion. Two half-sets of these images are then se- 
lected at random and subtracted from each other 
to produce one realization of noise. This technique 
is similar to the jackknifing b y scan commonly 



is similar to tne jacKkmnng b y scan c ommonly 
used for AzTEC data (see. e.g..ls"cott et al .112003 . 
2010t|Perera et alj|2008t Uustermann et al.ll2010l) . 



in that each jackknifed subset uses the scanning 
information of every available map. Use of this in- 
formation is especially important for our chopped 
data if we are to remove negative flux artifacts 
from the triple-beam PSF as well as positive flux. 
These "jackknifed noise maps" are more amenable 
to mass production, and are guaranteed to remove 
all contributions from a source however faint, so 
they are used in our pixel flux distribution (PFD) 
analysis (Q and in our Monte Carlo simulations 
to estimate numbers of spurious detections ( ^5.41) . 

The PFDs for S/N maps created with both 
jackknifed and shuffled noise exhibit random 
Gaussian noise to high precision (Figure [3]), with 
reduced chi-square for standard normal distribu- 
tion fits of 1.0 ± 0.2 and 1.2 ± 0.2, respectively 

3.3. Simulated maps 

Our simulated sky maps are constructed by 
populating noise maps with simulated sources. 
Careful construction of these maps is important 
for the fluctuation analysis described below ($4]), 
for which our method relics entirely on our abil- 
ity to authentically reproduce the signal from 
the MAMBO array so as to faithfully repro- 
duce the PFD. Thus, when adding sources into 
a noise map, we need to take into account the 
position-dependent negative sidelobes as well as 
the position-independent positive flux profile for 
each injected source. 

To handle the varying PS F properly, we take 
an approach similar to that of Greve et al. ( 20081 ) 
and model the changes in the PSF explicitly as a 
function of position. We use the MOPSIC pipeline 
script map_negres .mopsic, which will calculate 
the expected negative residual pattern on the sky 
in equatorial coordinates for a given set of obser- 
vations and an ideal, gridded, input source model. 
As an input we used an array of ideal Gaussian 
point source profiles, each with 11" FWHM, span- 
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ning the entire field and spaced as closely as pos- 
sible without having the sidelobes overlap. This 
minimum spacing is set by our larger chop throw 
(42" for all of our "best" data and the overwhelm- 
ing majority of our "full" dataset). The result is 
an array showing the full PSF near any location in 
the map (see Figured]). Because the ~ 84" spac- 
ing is less than the 300" extent of each individual 
map and the 120" separation between map point- 
ing centers, the PSF morphologies change slowly 
from one to the next. We thus generate an au- 
thentic point source response in a simulated map 
by using the closest available PSF relative to the 
position of a given injected source. 

4. P(D) analysis of the pixel flux distribu- 
tion 

We constrained the 1.2 mm number counts be- 
low our nominal sensitivity and confusion limits by 
performing a fluc tuation analy sis, also known as 
a P(D) analysis (ICondonlll974h . The P(D) anal- 
ysis has the advantage of using information from 
the entire PFD of the map (see Figure [5]) to con- 
strain the number counts, rather than using only 
those pixels above the source detection threshold 
(e.g., by counting bright sources). This distinction 
makes the P(D) analysis robust against the small 
number statistics of counting detections in the 
map. Additionally, the nature of the Monte Carlo 
simulation described below allows us to minimize 
uncertainties in flux boosting and completeness, as 
well as the effects of confusion and source blend- 
ing, because they are built into the simulation 
through the injection of model sources. These 
benefits have l ed Monte Carlo simu l ation PfD) 
analy ses (e.g., iMalonev et al.l l2005t IScott et al. 
20101 ). as well as Markov Chain Monte Carlo 
Metropoli s-Hastings (MCMCMH ) PfD) analy- 
ses ( e.g., Patanchon et al.l 120091 : IValiante et al 



20ld : Idenn et al.l l2010h . to be applied to both 



chopped and unchopped data at many wave- 
lengths. Our implementation of a Monte Carlo 
simulation PfD) analysis, which is best suited to 
handle our position-depende nt PSF, adopts the 



methods of lScott et all (|2010h 



The basic approach of our P(D) analysis is 
to parametrize the differential number counts 
and add a simulated map of sources obeying 
these number counts (along with their position- 



dependent negative sidelobes; see ^I3.3|) to a jack- 
knifed noise map (see §3.2|) . thereby creating a 
fully simulated MAMBO sky image. Because of 
its simple form, our initial parametrization is a 
single power law with normalization N^ m j y and 
index S, such that the differential number counts 
have the form 



dN _ 



4 mjy 
S 



(1) 



We adopt this form from lLaurent et al.l (|2005[ ). so 
as to minimize the degeneracy between the nor- 
malization and slope of the number counts at the 
flux density of our typical significant detection 
(~ 4mJy). Next, we optimally filter this fully 
simulated sky image and compare its PFD to that 
of the real data using the likelihood (see below) as 
a goodness-of-fit statistic. We then compute the 
average likelihood of the data for ten iterations 
of these model parameters, choose new parame- 
ters, and repeat the process. After filling parame- 
ter space with likelihood statistics, we identify the 
best-fit parameters as those giving the maximum 
likelihood. After the location of this peak in pa- 
rameter space is identified, we return and sample 
this one position ~ 10 3 times in order to constrain 
the absolute likelihood value enough to discrimi- 
nate between fits using different flux density cut- 
off values (see below) . 

The likelihood for each sky realization is calcu- 
lated as follows. Assuming the PFD's flux bins are 
uncorrelated, the probability of observing m pix- 
els in the ith flux bin given an expectation value 
of Xi is given by a Poisson distribution: 



A™* 



(2) 



Therefore, the natural logarithm of the probability 
P[ni\ of observing the full PFD {n^} for a model 
PFD {Xi} (the log- likelihood) is given by 

lnP[nJ = JZ ln (~ — i — ) = ^2n i lnX l -lnn i \-X l 

(3) 

We limit the comparison to bins in which rij > 10 
and use Stirling's approximation to write the sum 
as 

lnPM ~ ~ A i - n i ln (y) ( 4 ) 
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In reality our histogram bins are not uncorrelated, 
since our beam solid angle is ~ 10 x larger than 
the area of one pixel; thus, this expression will 
serve simply as a comparative metric for choos- 
ing a set of best-fitting parameters. The prop- 
erly calibrated error bars for this estimate can 
then be found via a Monte Carlo simulation us- 
ing synthetic images generated from the best-fit 
model. P[rii\ is therefore a function over the two- 
dimensional parameter space of (AT4 m j y ,5), within 
which the best-fitting model parameters are those 
that minimize — lnP[nj]. 

We apply the P(D) analysis to the region in 
the "best" map where the local RMS noise a < 
1.25 mJy beam -1 before filtering (see Figure [1}. 
This threshold was chosen to maximize the dis- 
criminating power of the simulation. If the noise 
threshold is very low, the region used for analysis 
has very high sensitivity, but there are fewer pixels 
available for comparison. If the noise threshold is 
too high (e.g., we use the full extent of the "best" 
map), too many regions with differing local noise 
properties are included, and the signal from the 
interior of the map is washed out. Our choice of 
threshold represents a compromise between these 
two limits, including as many pixels in the analy- 
sis as possible while keeping their noise properties 
as uniform as possible. 

Our initial expression of the PFD in terms 
of 15 bins between —3.5 and +4.0 mJy beam" 1 , 
chosen so that all flux bins had > 10 pixels, 
failed to constrain the model parameters with a 
unique maximum likelihood. Because the bright- 
est (and most model-constraining) pixels in the 
histogram are in the sparsely populated bins above 
4.0 mJy beam -1 , our simulations produced only 
a best-fit arc in parameter space. When we in- 
creased the sensitivity to the brightest pixels by 
adding an additional bin spanning from 4.0 — 
5.0 mJy beam -1 , wide enough to include > 10 pix- 
els, a well-defined global maximum appeared along 
the previously degenerate arc. 

In order to keep the models from diverging 
at low flux densities, we also imposed a faint- 
end cutoff in flux density, Scut, which we crudely 
treated as a third parameter in the P(D) analy- 
sis. By stepping through the values S cu t = 0.3, 
0.2, 0.1, 0.05, and 0.01 mJy, testing each with a 
full set of fitting parameters, we found the max- 
imum likelihood values for the power law to be 



41.7±0.4%,47.7±0.2%, 52.1 ±0.2%, 53.0±0.3%, 
and 49.7 ± 0.5%, respectively. The fits improved 
steadily with decreasing S cu t down to 0.05 mJy 
but then worsened at 0.01 mJy. Thus, the over- 
all best fitting parameters for the power law were 
A^mjy = 19.7±^deg- 2 mJy-\ 6 = 3.14±°;1|, 
and S cu t = 0.05 mJy (we quote marginalized 68% 
double-sided error bars). 

Our second num ber counts parametrization was 
a Schechter ( 19761 ) function of the form 



^=N> 
dS 4mJy 



4mJyV 
S J 



exp 



S — 4 mJy 



(5) 

Because the full four-dimensional parameter space 
of the Schechter function (N^ mJ , 5', S' exp , S cu t) 
is too large to probe with a blind grid-searching 
routine, we began by fixing S' cut equal to the solu- 
tion for the power-law model (0.05 mJy). We then 
alternated the P(D) analysis between varying the 
parameters (A^ mJy , S cxp ) and (A^ mJy , 5') until 
the solutions converged on the same values for all 
three parameters (convergence was achieved after 
three iterations). The initial seed guess for S' exp 
was motivated by naivel y scaling th e PfD ) solu- 
tion at 1.1mm, found bv lScott et al.l ( 2010l ) using 
AzTEC data in the GOODS-S field (S exPi lA mm = 
1.30 mJy), to 1.2mm (see H7.1[) . The results con- 
verged to the exponential scaling flux density of 
'S'exp = 1-05 mJy, which was then held fixed while 
we proceeded to make full searches over the pa- 
rameters (N' imJ , 5') while varying S' cut . 

The quality-of-fit for the Schechter function 
was also greatest for a flux density cutoff of 
S cut = 0.05 mJy. The maximum likelihood val- 
ues for S cut = 0.2, 0.1, 0.05, and 0.01 mJy 
were 23.8 ± 0.1%, 26.9 ± 0.1%, 28.3 ± 0.2%, and 
26.7 ± 0.4%, respectively. The final set of best- 
fitting parameters for the Schechter function were 



14.5^7 deg -2 mJy- 



1.86 



+0.20 
-0.23' 



/V' 

JV 4mJy 

S' exp = 1.05 mJy, and S cut = 0.05 mJy. 

Figure [6] shows full 68% and 95% confidence 
regions around the parameters of maximum likeli- 
hood for both the power law and Schechter func- 
tion parametrizations. The uncertainty contours 
were ge nerated via Mont e Carlo sampling (used, 



e.g., m IScott et al.l l2010h by taking the best-fit 



model parameters and using them to generate ad- 
ditional simulated sky realizations. A P(D) anal- 
ysis was then carried out on each of these rcaliza- 
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tions to recover a set of new (scattered) best-fit 
model parameters. This process was performed 
^100 times with the same model inputs in order 
to generate a likelihood density map around the 
best-fit model parameters. 

As an additional constraint on our model fitting 
and a means of choosing between paramctriza- 
tions, we also require that the number counts 
model obey the constraint of the 1.2 mm CIB (see 
> j7.5l) , which is shown as the shaded region in Fig- 
ure [6] It is evident in Figure [6] that although the 
power law parametrization can fit our observations 
well, it significantly overpredicts the 1.2 mm CIB. 
In contrast, the Schechter function parametriza- 
tion is in excellent agreement with the constraint 
of the CIB. We therefore adopt the Schechter func- 
tion parametrization as our fiducial model (e.g., 
Figure [5]) for the remainder of the paper. At the 
highest flux densities the Schechter function and 
power-law models nominally predict very differ- 
ent behavior; however, our fluctuation analysis is 
not sensitive to the number counts at flux den- 
sities higher than those of our brightest detected 
sources. 

The fact that both the power-law and Schechter 
function models of the differential number counts 
fit best when S cut = 0.05 mJy suggests that the 
1.2 mm number counts do not keep rising far be- 
yond 0.05 mJy; formally, they may begin to fall 
between 0.05 and 0.01 mJy, or may already be 
decreasing by 0.05 mJy. The former case is in 
agreement with recent survey s of lensing clus- 



ters (e.g., iKnudsen et all 120081 ) that show SMC 



number counts increase at least as far down as 
~ 0.1 mJy at 850 fim. At 1.2 mm, this corresponds 
to Si. 2 mm — 0.03 — 0.04 mJy using the submil- 
limeter spectral indices determined from match- 
ing detections at 850 fxm, 1.1mm, and 1.2 mm in 
the GOODS-N and COSMOS fields (Greve et al. 
2008; Chapin et al. 2009; Austcrmann et al. 2010; 
se e ^7.11) This result is in contrast to the analysis 
of lScott et a"fl(|2010l) . who found that the choice of 
S cu t did not affect their results. The discrepancy 
may be due to the fact that their lower resolu- 
tion (28" HPBW) reduces the effective depth that 
can be reached before sources begin crowding in 
the beam, thereby reducing the sensitivity of the 
P(D) technique. 



5. Analysis of bright sources 
5.1. Source extraction 

We extracted sources from our "best" and "full" 
maps in a three-step process. First, we minimized 
the chi-squarc statistic for a two-dimensional 
Gaussian profile with an 11" FWHM fit at each 
pixel center. This minimization was achieved 
quic kly by using a matche d filter convolution (see, 



e.g., Serjeant et al. 20031) . Given a signal image 



Sij, an image of the local RMS noise tjjj, and a 
smaller array describing the telescope's PSF P xy , 
the chi-square statistic for a source with flux F 
located at pixel can be written 



X 2 (F\iJ) 



( Si—x,j—y FP xy 

xy 



&i-x,j-y 



(6) 



We ignored the position-dependent negative side- 
lobes when we applied the matched filter and used 
only the central Gaussian profile, since the nega- 
tive residual flux reaches only < 5% of the peak 
positive intensity (see Figure 0]) in the map inte- 
rior. Additionally, our significant detections are 
on average farther away from each other than the 
largest chop throw used during the observations 
(48"), so their effect on our flux measurements will 
be less than 5% of our strongest sources' flux den- 
sities (i.e., < 0.25 mJy). 

By finding the minimum of x 2 (P\h j) as a func- 
tion of F a nd determining the associated uncer- 
tainty AF (jSerieant et al 1 120031) . we obtained 



dF 







F 

Xf 



^2xy & 



■x,J-y 



Wi- 



x,J 



-yPxy 



x,J 



p2 
-y 1 xy 



(7) 

as the S/N of each pixel, in terms of the weight 
map W defined in £13.11 Next, we located the 
centroid of each source to sub-pixel precision by 
fitting the PSF to the region in the original sig- 
nal map at the location of each significant peak in 
the S/N map, allowing the position of the Gaus- 
sian to vary. Figure [7] shows the uncertainty in 
this best-fit centroid position, derived via Monte 
Carlo simulations. For the typical flux densities 
of our significant detections, the average offset be- 
tween injected and recovered centroids is 1" — 3". 
Finally, we computed the best-fit flux density by 
taking the matched-filter weighted average of the 
flux map, this time with the PSF kernel shifted 
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by interpolation to the more precise location of 
the source centroid. After the flux density of each 
source was recorded, the source was removed from 
the map by subtracting the flux-scaled PSF from 
the source location before we searched for the next 
most significant detection. 

By propagating the uncertainty in the sig- 
nal image thr o ugh t he ^-minimization process, 



Serjeant et al.l (|2003l ) have shown that the un- 
certainty in the resulting best-fit flux density at 
position is 



AF(i,j) = 



1 



p2 

-y xy 



(8) 



We find that this expression consistently under- 
estimates the uncertainty in our map. The rea- 



son is simply that the derivation bv lSerieant et al 



implicitly assumes that the Gaussian noise 
in the flux image is spatially uncorrelated. In our 
images, the noise is correlated on a length exactly 
matching the FWHM of the telescope PSF, and 
pure noise fluctuations can be amplified along with 
the real point sources. In order to correct for this 
underestimate we (1) produce an optimally filtered 
image and weight map, (2) use this filtered image 
and weight map to make a S/N ratio map, and 
(3) rescale the filtered weight map such that the 
standard deviation of the S/N map equals unity. 
This empirical calibration corrects for the effects 
of correlated noise and produces a map with accu- 
rate post-filtered flux density uncertainties. The 
accuracy of the method is evidenced by its match- 
ing the predicted number of positive excursions in 
a correlated Gaussian field as a function of S/N 
level (see HSU}. 

The quoted Ict errors (see columns S®° st and 
S^ ul1 in Table [5]) for each detection correspond to 
the rescaled version of AF evaluated using the 
shuffled noise map at the location of the source 
(see §3.2j) . This noise map describes the Gaus- 
sian noise of the observations more faithfully than 
the original signal map, which overestimates the 
noise by ~ 5% due to the positive and negative 
sidclobcs from bright sources. Because our noise 
is well above the estimated confusion limit ( £I5.3|) . 
Gaussian random fluctuations are the dominant 
source of uncertainty in our measurements. 



5.2. Comparison of results for best and 
full maps 

As discussed in control software problems 
during our first two semesters of observations un- 
dermined our confidence in the reliability of the 
resulting maps. To assess whether the "full" 
map could be trusted for bright source extrac- 
tion, we performed two comparisons between our 
"full" data and observations with pristine calibra- 
tion. For the first comparison, we carried out the 
source extraction steps described in Section 15. fl 
for both the "best" and the "full" datasets and 
compared the properties of the sources recovered 
from each. Specifically, we began by choosing the 
eight sources with S/N > 4.5cr detections in our 
"best" map: above this threshold, we expect to 
see fewer than one spurious detection ( §5.4|) . All 
eight of these sources are recovered with > 5. Oct 
significance in the "full" map. Figure [5] shows the 
locations of these sources in the field, and the lo- 
cations of the additional > 5. Oct detections in the 
"full" map. Each of the eight sources increases 
in significance between the "best" and the "full" 
maps. Additionally, all but one of the "full" map's 
17 > 5.0ct sources are identified in the "best" data 
at lower significance. One source kept the same 
significance because it lies in the northeast corner 
of the field, where observations in the first two 
semesters contribute little additional sensitivity. 
Further, for these 17 sources, the ratio of the flux 
densities in the "best" and "full" maps is consis- 
tent with unity (Figure [9]). 

Next, we compare d our "full" map to MAMBO 
on-off photometry bv lFiolet et al. ( 2009 ) of Spitzer- 
selected high-redshift starburst candidates in the 
LHN. We tabulated the map flux densities at the 
positions of the 13 galaxies in their sample that lie 
within our map's footprint (two of these turn up 
as significant detections in our "full" map; we use 
these sources' non-deboosted flux densities here) 
an d compared them to the flux densities reported 
bv iFiolet et~aH (|2009h . We found that the flux 
densities from the two significant detections as 
well as those from 10 of the 11 non-detections are 
consistent to within Ict (see Figure [TOJ one source 
is only consistent to within ~ 1.5ct). 

These two successful consistency checks lead 
us to conclude that the errors in our first two 
semesters' data are not at a level that compro- 
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mises point source detection, at least for high- 
significance sources. We have therefore proceeded 
to define our source catalog on the basis of the 
"full" map. Since the fluctuation analysis decribed 
in fJ3] relies on the authentic reproduction of the 
field's noise properties and low-S/N fluctuations, 
we have restricted this analysis to the "best" data 
only. 

5.3. Confusion 

Random Gaussian noise is the uncertainty in 
the total flux density inside any single beam on 
the sky due to random fluctuations, while confu- 
sion noise is an additional uncertainty in the flux 
density of a single source due to the contributions 
of faint sources within that beam. The "confusion 
limit" is defined as the flux density threshold at 
which confusion noise significantly affects the mea- 
sured flux density of a source, and is commonly 
taken to be the flux density above which the inte- 
grated number coun ts of all brighte r sources reach 
~ 0.033 per beam ( Condon! 1974 ). In our map, 
this rule gives ~ 0.9 mJy using 9b = 15.6" (in the 
smoothed -H> filtered version of our "full" map) 
and assuming our best- fit number counts (Q. 

We have also made a direct estimate of the con- 
fusion noise by comparing the noise in the central 
regions of the filtered "best" map to the same re- 
gion in a series of filtered jackknifed noise maps. 
Since the jackknifed noise maps remove confused 
as well as bright sources, the increase in average 
RMS noise in this relatively uniform region in- 
dicates our map contains confusion noise at the 
level of o~c — 0.24 mJy beam" 1 . As a consistency 
check, we have also estimated the confusion noise 
by generating simulated maps with source popu- 
lations following our best-fit number counts from 
0.05 mJy up to the confusion limit of 0.9 mJy. Due 
to the central limit theorem, these faint and con- 
fused maps with zero mean have roughly Gaus- 
sian PFDs and provide approximations of the con- 
fusion noise, assuming our model of the number 
counts. The standard deviation in these maps 
is 0.21 mJy beam -1 , in agreement with the mea- 
sured confusion noise within the uncertainties of 
the number counts model. We therefore adopt the 
measured value of ctc — 0.24 mJy as our estimate 
of the confusion noise. The average uncertainty 
in the flux density of a source in our catalog is 
0.62 mJy, indicating that confusion does not dom- 



inate our noise budget. 

5.4. Spurious sources 

We estimated the number of spurious detections 
as a function of S /N by running our source extrac- 
tion algorithm on various noise maps (see §3.2j) . 
We tested jackknifed noise maps, shuffled noise 
maps, and simple Gaussian random numbers. The 
Gaussian random numbers had a spatially varying 
standard deviation matched to the weight map of 
the observations. Figure [Til shows the mean total 
numbers of spurious detections found in 10 3 jack- 
knifed and Gaussian number noise maps, and in 
10 2 shuffled noise maps, as a function of S/N. All 
three styles of noise map are consistent with each 
other in their ability to produce spurious detec- 
tions with S/N > 3.0. This result confirms that for 
the purposes of extracting high-significance detec- 
tions, the shuffled noise maps are just as "source- 
free" as the jackknifed noise maps. Additionally, 
both are consistent with a Gaussian distribution 
down to 3. Oct (and likely consistent with Gaus- 
sian noise at all S/N, as implied by the PFD his- 
tograms in @. The over-plotted curve in Figure 
[TT] shows the expected number of excursions above 
a given S/N level in any isotropic and homoge- 
nous Gaussian random field, derived (and thus 
only formall y valid) for h igh excursions (see, e.g., 
Chapter 6 oflAdl er Ill98lh . The agreement at high 
S/N indicates that our noise maps and source ex- 
traction algorithm are well-behaved. For both the 
"best" and "full" maps, we expect 0.8 (5.4) spuri- 
ous sources with S/N > 4.5 (4.0). 

A source at high risk of being a spurious de- 
tection can be identified by calculating the to- 
tal probability that th e deboosted flux density 
is < OmJy (see , e.g., Austermann et all 120101 



IScott et a"i1l2010h . hereafter referred to as P(S < 
0). U sing the threshold of P( S < 0) > 0.10 
used bv l Austermann et al. ( 2010f ). we identify only 
one (ID # 33) high-risk spurious detection in our 
S/N > 4.0 sample (see Table 

5.5. Completeness 

We estimate the completeness in our data us- 
ing the Monte Carlo method of searching for in- 
jected sources of varying flux density. The inho- 
mogeneous noise in our map means that sources of 
identical flux density have different probabilities of 
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being detected in different locations. We account 
for this effect statistically by performing the com- 
pleteness simulation assuming sources have equal 
probability of being located anywhere in the map 
during the injection process. Although the high- 
redshift star-forming galaxy population that our 
observations trace is likely to exhibit clustering, 
with the brightest galaxy mergers occurring in 
the most massiv e dark matter halos (see, e.g., 
Weifi et al.ll2009t ). the large Az interval to which 
millimeter selection is sensitive tends to weaken 
angular clustering signatures (see however i j5.8l) . 
We inject model sources with known flux den- 
sity into our original signal map one at a time at 
random positions and search for them using the 
same source extraction algorithm used to create 
our source list. If an artificial source is recovered 
with S/N > 4.0 within 11" of the injected loca- 
tion, it is considered detected. The injection pro- 
cess was repeated 10 3 times for each flux density 
in a logarithmic grid from 1.0 mjy to 10.0 mjy; 
the average recovery percentage is shown in Fig- 
ure Q21 Our map is 80% complete at 3.7 mjy and 
50% complete at 2.6 mjy. We also tabulated the 
angular separations between the injected and re- 
covered source positions to characterize the un- 
certainties in the positions of our actual signif- 
icant detections (ignoring telescope pointing er- 
rors). For Si. 2 mm > 2mJy, the average position 
error (A0) < 3" (see Figure EJ}. 

5.6. Flux boosting 

To correct the measured flux densities of our 
detections for the effect of "flux boo sting," we use 



the B ayesian technique described in iCoppin et al 
( 20051 ). which because of its versatility in han- 
dling both chopped and unchopped data has 
been adapted for use at many wavel e ngths (e.g., 
Coppm et aDl2006t iGreve et all 12005 IScott et al 



P(S,a\S ) is the likelihood, and P(S,a) is the 
prior measurement distribution. P{S, a) is in- 
dependent of So, so only acts to normalize the 
expression such that J P(So\S,a)dSo = 1; here- 
after, it will be ignored. We have shown that the 
uncertainty in our map is dominated by Gaussian 
random noise, so P(S, cr\So) takes the form 



P(S, <t\Sq) oc exp 



(S - So) 2 
2a 2 



(10) 



To estimate the prior flux distribution in the map 
(P(S )), we assembled a PFD containing the pix- 
els from 10 4 noise-free random sky realizations 
that used the best-fit number count parameters 
from our P(D) analysis. The peak value and 68% 
double-sided confidence intervals of the resulting 
posterior probability function (P(So\S,a)) were 
found numerically for each measured flux density 
and uncertainty. Figure [T3] shows four examples of 
the deboosting process, each in a different regime 
of source S/N. If the S/N is too low, and the in- 
tegration of the confidence intervals does not con- 
verge, we instead use an analytic formula to es- 
timate the deboosted flux density. For this, we 
generalize the formalism of lHogg &: Turner dl998h 
to a Schechter function, and locate the maximum 
of the posterior flux distribution: 



P(Sq\S, a) kS s exp 



Sp 
SL 



(S-Sp 
2a 2 



(11) 

where 5' and S' exp are the power-law slope and ex- 
ponential scale factor of the Schechter function, 
respectively. By solving for So when the deriva- 
tive of the above expression vanishes, we find the 
highest posterior probability to be achieved for 



2008 1 ; IPerera et al.l2008HScott et al.l2O10HAustermann et 
20101). Using Bayes's Theorem and the prior in- 



formation of the number counts functional form 
found from our P(D) analysis O, the probabil- 
ity that a source has true flux density So given a 
measurement S with uncertainty a is equal to 

P(S,a\So)P(S ) 

PS S,cr = — — 9 

P{S,a) 

where P(S,<j\Sq) is the posterior probability, 
P(Sq) is the prior flux density distribution, 



al. 

Struc 



'a 2 -SS, 



exp ) 



-4 5>S>l p a 2 



2S' 



(12) 

To ensure that our adaption of the Bayesian 
method of flux deboosting returns a properly cor- 
rected estimate of the true number counts, we 
performed a Monte Carlo simulation to directly 
calculate the observed number counts of random 
sky realizations populated with source distribu- 
tions following our best-fit number counts (see, 
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e.g., ICoppin et aLll2006f) . Figure [14] shows the re- 
sults of this consistency check. This simulation 
demonstrates that the Baycsian method of flux de- 
boosting performs well in recovering the original 
injected number counts. The residual scatter of 
the average recovered number counts around the 
average input model in Figure[T4]demonstrates the 
level of systematic error in the algorithm, which is 
significantly smaller than the statistical error of 
our differential number counts estimate (see Fig- 
ure [15]). 

5.7. Direct calculation of number counts 

While our catalog of detections includes all 
sources with S/N > 4.0, we use only detections 
with S/N > 4.5 for our direct calculation of the 
number counts because above this threshold, we 
expect to detect less than one spurious source (see 
Figure QT]). Table [3] presents integral and differen- 
tial number counts after correction for complete- 
ness and flux boosting. Figure [15] shows our di- 
rectly calculated number counts, along with the 
95% confidence regions for the best-fit power law 
and Schcchtcr function models of the differen- 
tial number counts found from the P(D) analysis. 
These two independent methods of estimating the 
number counts are in agreement with each other. 
This consistency is encouraging because the P(D) 
analysis and the direct estimate of number counts 
depend on the faint and bright pixel values in dif- 
ferent ways. 

5.8. Clustering 

The group of sources in the southeast corner of 
our field, as well as the large void in the center, 
prompted us to perform a clustering analysis to de- 
termine whether or not the distribution of sources 
in our map is statistically clust ered or not. To 
perfor m the analysis, we used the lLandv fe Szalav 
(|1993l ) correlation function estimator: 

DD - 2DR + RR 
w ( d ) = ( 13 ) 



with variance 



(l + w(fl)) s 
RR 



(14) 



( Gawiser et al ] l2006h . In the equations above, 
DD, RR, and DR represent the normalized num- 
bers of unique galaxy-galaxy, random-random, 



and galaxy-random pairs with angular separations 
9 ± d6/2. This estimator is used freq uently in ex- 
tragalactic deep field analyses (e.g.. iBorvs et al 



2003tlScott. Dunlop. fe Serieant|[200a IWeifi et al 



2009), and has been shown to have nearly Pois 



son variance and zero bias (|Landv fe Szalavill993l ). 
We take into account the geometric boundary of 
the map and the variation in detectability with 
position by generating the random locations with 
the same Monte Carlo algorithm used for the P(D) 
analysis. We inject ensembles of sources following 
our best fitting number counts into a noise map at 
random locations and use the positions of sources 
detected with S/N > 4.0 as our random coordi- 
nates. This Monte Carlo technique is important 
for ensuring that we do not misinterpret depth 
variation in the map as a clustering signal. 

To confirm that this technique is unbiased, we 
also performed the full clustering analysis on only 
random positions to check that we recovered a flat 
w{9) = response (see Figure [16]). For small sep- 
arations (< 2'), however, it turns out that w(9) 
does not return zero in our data: depending on 
the position within the map, the negative sidelobcs 
can suppress the flux densities of nearby sources 
enough to lower their S/N ratios below the detec- 
tion threshold. This effect begins to have an effect 
at ~ 2 x the chop throw (of which the maximum 
used in any semester was 48"), and has a strong 
effect at separations < lx chop throw. Because 
this effect suppresses the detection of RR and DD 
pairs but not DR pairs, the zero-clustering base- 
line for chopped data like ours is less than zero 
at these small angles. In order to assess the clus- 
tering in the map while taking this bias into ac- 
count, we measure the effective clustering relative 
to the zero-clustering baseline for these separa- 
tions (< 2'). 

The result of our Monte Carlo clustering anal- 
ysis is shown in Figure 1161 We find a small 
clustering signal when using all detections with 
S/N > 4.0 that agrees reasonably well with 
the angular correlation f unctio n measured by 
IScott. Dunlop. fe Serjeant! (|2006h . who combined 
many different SCUBA 850 /im blank field maps, 
and that shows stronger clustering (albeit at lower 
S/N) than t he cor relation function measured by 



Weifi et all (l2009h a t 870 /xm in the ECDF-S. 



Williams et al.l (|201ll ) have analyzed the cluster- 
ing of 3.0-3.5cr 1.1mm detections in a 0.72 deg 2 
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map of the COSMOS field with ASTE/AzTEC, 
concluding that it is difficult to recover reliable 
clustering parameters for SMGs from maps whose 
angular resolution and total area are limited. This 
result argues for caution in interpreting our clus- 
tering analysis, although we do benefit to an ex- 
tent from MAMBO's relatively high angular res- 
olution. An interesting feature in our correlation 
function is the spike near 9 ~ 4'. This signal is due 
to the rich group of sources in the southeastern 
corner of the map, all at typical relative spacings 
of a few arcminutes from e ach other (se e also £ 17.41) . 
It may be noteworthy that lWeiB etaD (|2009h find 
a ~ 2.4ct spike above their best f itting model of 
l o(9) a t a scale of ~ 5', near where Williams et al.l 



(|201u) also detect a slight positive excess in oj(6). 
When performing the analysis on only our (27) 
most significant sources with S/N > 4.5, we find 
no significant clustering signal. 
Following the analysis of, e.g 



Borvs et al 



(120031) for the SCUBA "Supermap ," we also use 
the method of lScott fe Tout! (|1989|) to analyze the 
cumulative distribution of nearest neighbors to 
test whether our galaxy positions are consistent 
with being drawn from a random distribution (see 
Figure [TT]) . Because the nearest neighbor anal- 
ysis is sensitive to the total number of positions 
used, we use the 41 most significant detections 
in each Monte Carlo realization, instead of all of 
those detections with S/N > 4.0 as in the cor- 
relation function analysis. Because the number 
counts rise quickly, the S/N of the least signifi- 
cant discovered source varies, but is always close 
to 4.0. A Kohnogorov-Smirnov test rules out the 
null hypothesis that our significant detections are 
drawn from a random position distribution at the 
95% confidence level, implying that the source 
locations in the map (e.g., defining the southeast- 
ern clump and the central void) are not arranged 
randomly. 

6. Counterpart identification 

We have calculated th e corrected probabili ty of 
chance associations (P; iDownes et al. 19861 ) be- 
tween our MAMBO detections and possible coun- 
terparts at the other wavelengths at which the 
LHN has been observed (see Table 0). The P 
statistic is defined by 



P = l-e~ E (15) 

for E = P*[l + ln(P c /P*)], P* = Trr 2 N(> S), 
and P c = irr^N c , in terms of the brightness of the 
counterpart 5 1 , the source separation r, the search 
radius r s , the number density of sources brighter 
than S N(> 5), and the number density of sources 
at the sensitivity limit N c . Based on the results 
of the position error analysis (see Figure [7]), we 
chose a counterpart search radius of 8". Because 
positional uncertainty a cx FWHM x SNR -1 , our 
11" beam is the dominant source of error, and we 
ignore the positional uncertainties at other wave- 
lengths. We consider P < 0.01 to define a robust 
counterpart, 0.01 < P < 0.05 a likely counterpart, 
and P > 0.05 an unlikely association. 

6.1. 20 cm counterparts 

We used the NRAO VLA Sky Sur vey (NVSS) 
and th e de ep SWIRE number counts of lCondon et al 
(|l998l ) and lOwen fc Morrisonl (|2008h . respectively, 
in the calculation of P to assess the significance 
of 20 cm counterparts. The 20 cm VLA point- 
ing of the LHN has a central RMS sensitivity of 
2.7 /iJy, rising to ~ 4-5 /iJy near the edges of 
our MAMBO map. When we compare our 41 
S/N > 4.0 detecti o ns to the 5cr 20 cm catalog of 
Owen fe Morrisonl (|2008h . 44% (18) have robust 
counterparts, and 41% (17) have likely counter- 
parts. We have also reexamined the 20 cm map 
in the vicinity of the remaining MAMBO sources 
and have identified one additional robust coun- 
terpart (ID #9), two likely counterparts (ID #28 
and # 36), and one unlikely counterpart (ID # 20) 
at the 4-5cr level. We also deblended one likely 
counterpart into one robust and one unlikely coun- 
terpart (ID 17). After including these additional 
sources, 49% (20) of our MAMBO detections have 
robust counterparts, 44% (18) have likely coun- 
terparts, and 7% (3) have unlikely or no detected 
counterparts. We performed a Monte Carlo sim- 
ulation to test the reliability of our P values and 
found that 4.9 ± 0.2% of randomly chosen posi- 
tions within our MAMBO field have a likely radio 
counterpart (P < 0.05) within 8", confirming the 
validity of the high number of robust associa- 
tions. We expect ~ 5 spurious detections above 
S/N > 4.0; thus, the handful of sources with un- 
likely or no radio counterparts may be spurious 
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detections if they do not lie at a very high red- 
shift. One MAMBO source (ID #32) with an un- 
likely (fbocm = 0.056) radio counterpart also has 
a likely (P250um = 0.016) Herschel counterpart 
( Magdis et al. 201fj() . arguing against its being a 
spurious detection. 

6.2. 50 cm radio counterparts 

We have extracted 50 cm flux densities from 



the GMRT map (jOwen et all l2011al) with the 



same technique used at 20 cm and 50 cm (see 
Owen fc Morrisonll2008t lOwen et al.ll2009bl ). The 
uncertainties listed in Table [5] reflect the local 
RMS noise in the image and do not include a ~ 3% 
calibration error or a spatially varying GMRT 
pointing error. Two of the 50 cm detections are 
heavily blended with bright neighbors, so for these 
counterparts we report only tentative fluxes. Of 
the 40 tabulated 20 cm counterparts (including the 
two with P > 0.05), all 40 have 50 cm counter- 
parts. The one 20 cm non-detection (within 8") is 
also a 50 cm non-detection. 

6.3. 90 cm radio counterparts 

To search for 90 c m counterpar t s, we u sed the 
90 cm radio catalog of Owen et al. ( 2009bl ) , which 
has an RMS sensitivity of 10 /j,Jy. Of our 41 
MAMBO sources, 24% (10) have 90 cm counter- 
parts. Each of the ten 90 cm counterparts is also 
detected at 50 and 20cm with Poo cm < 0.05. 

6.4. 24 /im counterparts 

In addition to SWIRE 24 yum observations of 
the LHN (3cr depth of 209 /iJy), there exist deeper 
Spitzer/MIPS data with a 3<r depth of 18 /J,3y 
( Owen et al ] l2011bh . We searched for 24 /im 
counterparts to our MAMBO detections in this 
deeper MIPS image. To calculate P statistics 
for 24 nm count e rpart s, we used the counts of 
Bethermin et al. ( 2010h . A Monte Carlo simula- 
tion of the 24 /im P-statistic finds that 4.7 ± 0.3% 
of random positions yield a counterpart with 
P < 0.05. Within our sample of 41 MAMBO 
sources, 20% (8) have robust 24 ^m counterparts, 
and 29% (12) have likely counterparts. 

6.5. X-ray counterparts 

Only one source (MM J104522. 8+585558 = 
ID # 26) has a likely X-ray counterpart (CX- 



OSWJ104523.6+585601: fWilkes et al.l l2009h. The 
X-ray source is at a distance of 7.2" and has a 
broad band (0.3-8. keV) flux of (2.5 ± 1.1) x 



10 15 erg cm 2 s 



(jPolletta et al.ll2006l) . By us 



- 1 ^ — 9 

° erg cm s 



ing the 2.5-7 keV flux of 1.58 x 10~ 
together with the Chandra . /SWIRE counts from 
2-8 keV (| Wilkes et al.ll2009h we can set an upper 
limit on the probability of chance association of 
P < 0.02. 



7. Discussion 

7.1. Number counts vs. 
fields 



previous deep 



Pre vious deep su r veys a t 1 . 2 mm using M AMB O 



(e.g., iGreve et al.l 12004 iBertoldi et al.l 120071) 



have returned directly calculated 1.2 mm num- 
ber counts in the Lockman Hole East (LHE), 
ELAIS-N2, and COSMOS fields. The parame- 
ters of these surveys are listed in Table 0] we 
compare their results to our directly calculated 
counts, as well as to our best-fit P(D) models, in 
Figure [15] We find that our power-law slope is 
consistent with their results, but our results have 
a lower overall normalization. This difference in 
normalization might be due to the different meth- 
ods used in the number counts calculations. We 
have used the B ayesian metho d of flu x deboost- 
ing presented in ICoppin et al.1 (|2005h . and only 
include our most significa nt detections in the cal- 
culation. The analys es of Greve et al. ( 2004f ) and 
Bertoldi et al.l (|2007| ) use the method of injecting 
sources into noise maps to determine their flux 
deboosting correction, and include sources with 
lower S/N in their number counts calculation. In 
principle, any S/N cutoff would be acceptable for 
the latter calculation as long as the completeness 
correction uses the same threshold; however, lower 
S/N thresholds will lead to more spurious detec- 
tions. Both of these effects could be contributing 
to their higher normalization. However, consider- 
ing the relatively large error bars on all the mea- 
surements and the internal variation among the 
ELAIS-N2, LHE, and COSMOS datasets them- 
selves, the results are still nearly consistent. 

Because a single power-law parametrization of 
the number counts is commonly used to compare 
the results of deep surveys, we begin by noting 
that our best-fit power law index (S = 3.141" Jg) 
is consistent with the results of surveys at other 
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wavelengths that fit their number c ounts using a 



the GOODS-N field is S, 



simila r (single power-law) model. ICoppin et al 
( 20061 ) find 850 (im power-law indices of 8 = 
2.9 ± 0.2 and 8 = 3.0 ± 0.3 in the LHE and 
Subaru/XMM-Newton deep fields , respectively. 



(jPerera et al.1 12008c iChapin et al 



850 ,,m / S-\ 



mm ^ 1.8 - 2.0 

20091) . When 



comparing SCUB A and MAMBO detections, 
Greve et al. | (l2008h find Stsso^/Sx. 2mm — 

2.5. By 



Using Bolocam data at 1.1mm, lLaurent et al 
((2005J) estimate a power- law index of 8 = 3.16 
from directly c alculated counts in th e Lockman 
Hole. However, iMalonev et al.l (|2005l ) performed 
a P(D) analysis on the same 1.1 mm Bolocam data 
and find 8 = 2.7t„ 15- Alth ough it is well within 
the 1 a uncertainties of the Laurent et al.l ( 2005 ) 
result, the latter slope differs from ours by > 2<r. 
In this case, differences in the methods of our 
PfD) analyses might be the differentiating factor. 
Malonev et al. I (120051) also used chopped observa- 
tions in their analysis, for example, but ignored 
the effects of chopping on the PFD. It is possible 
that by not including the negative residual flux 
in their P(D) analysis, they required many fewer 
faint sources to match the pixel distribution of 
the real data (and therefore derived a shallower 
power law slope). However, a value of 8 ~ 2.7 
is also the bes t -fit si ngle power-law slope found 
by IScott et al.l (|2010h for their P(D) analysis of 
unchopped 1.1mm AzTEC data in the GOODS-S 
field. 

In Figure [15] (see also Table Q} , we also com- 
pare our results to those for deep field ob- 
servations at 1.1mm by AzTEC of the COS 



MPS (|Scott et al.l l 2008h. GOODS-N (IPerera et al 

2010|), 



2008 ). GOODS-S jlScott et al.l l2010l). SHADES 



coadding the MAMBO and AzTEC observations 
in the GOODS-N field into a map at an effective 
wavelength of A = 1.16 mm, iPenner et al.l (|201lh 
find an average value of Si.i6mm/ Si.imm ~ 0.88 
and Si.i6 mm /5i.2mm ~ 1.14. All of these results 
are consistent with a single modified blackbody 
spectrum, for /3 = 1.5 and Td = 30 K, observed 
at z ~ 2.5. Therefore, we adopt this fiducial 
galaxy model when comparing fluxes at differ- 
ent wavelengths and use 5850 ( um/<S'i.2mm = 2.3, 

5 l 870pm/5'l.2mm = 2.2, and Sl.lmm/Si^mm = 1.2. 

Our directly calculated counts are in excellent 
agreement with the rescaled results of the AzTEC 
surveys. Additionally, our prediction for the shape 
of the number counts below our sensitivity thresh- 
old, afforded by our P(D) analysis, agrees well 
with the deepest AzTEC number counts and is 
even in r ough agreement with the deepest SMG 
counts bv iKnudsen et al. ( 20081 ). 

Figure [15] also compares our results to vari- 
ous number count predictions derived from back- 
ward evolution models that incorporate multi- 
waveband observations of number counts and red- 
shift distributions, as well as limits imposed by 
the CIB light. We have restricted this compar- 
ison to models th at offer predictions at wave- 
lengths o f 1.2 mm (Bethermin et al. 2011) or at 



1 . 1 mm d Valiante et al 



(jAustermann et al.l 2010h. and AKA RI, SSA-N2, 120091 : iMarsden et al.ll2011l ). towhich we can apply 



20091 iRowan- Robinson] 



and SXDF (lHatsukade et al.l |2011|) deep fields 



(For c larit y, the observa t ions o f lHatsukade et al 



(|201ll ) and IPerera et al.l (|2008f ) are not shown in 
the plot because their data points lie within the 
scatter of the other AzTEC observations.) We 
also compare our results to the extremely deep 
SM G counts measured in lensing fields at 850 /im 
by IKnudsen et al. ( 2008 ) as w ell as the recent 
wide map bv IWeifi et all (|2009h using LABOCA 
at 870 fim in the ECDF-S. In order to compare 
our number counts directly to the results of these 
surveys at other wavelengths, we rescale their 
flux densities. Our choice of rescaling factor is 
based on the direct comparisons between Sssn/xm, 
Si.imm, and Si. 2mm for galaxies in the GOODS-N 
field. The average flux density ratio for sources 
with robust SCUBA and AzTEC detections in 



the rescaling described above. Although flux scal- 
ing will generally not provide a precise represen- 
tation of a model's predictions at 1.2 mm, the ex- 
trapolation from 1.1mm to 1.2 mm is fairly mod- 
est. At flux densities equal to or less than those 
of our significant detections, we find that our ob- 
servations are generally co nsistent with al l mode l 



predictions, although the Valiante et al. ( 2009^ 1 



model slightly overpredicts our P{D) curve near 
1 mJy. At the high flux density limit, all models 
uniformly ovcrpredict the counts from our best- 
fitting Schcchtcr function model, while remaining 
consistent with the predictions from our power- 
law result (which is only marginally compatible 
with the CIB; see Figure [5]). However, we cannot 
draw any conclusions from this apparent discrep- 
ancy, as our P(D) analysis cannot constrain the 
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differential counts at flux densities greater than 
those of our brightest detections. 

7.2. Fractional counterpart identification 

Here we investigate the question of whether our 
radio counterpart identification rate (-Ri. 4 ghz) 
of ~ 93±7% (38/41) in the LHN is intrin- 
sically greater than is seen in other surveys, 
or if it is simply a function of the increased 
20 cm sensitivity in this field. We compare 
our identification rate to those found for pre- 
vious deep surveys at 850 /im, 870 /im, 1.1mm, 
and 1.2 mm. Table [5] lists recent millimeter 
and submillimeter deep field surveys from the 
GOODS- N, LHE, SXDF, COSMOS , and ECDF- 
S fields feorvs et all 120031 |2004 llvison et al 



Perera et al.ll2008 
20091 ; iBiggs et al ' 



20071: iBertoldi et alJl2007t ISchinnerer et all 12007 



Chapin et al. 2009; Wei fi et al 



2011ft . along with their 20cm 



radio counterpart identification rates and 20 cm 
map sensitivities. Because the surveys have dif- 
ferent definitions of "significant" (sub)millimctcr 
detections, different data reduction techniques, 
and different standards for radio counterpart as- 
sociations, we marginalize over all of these vari- 
ables by looking at the average radio counterpart 
identification rate, and the average 20 cm map 
sensitivity. Using the surveys listed in Table we 
find (cti.4ghz) ^ 7.2 /iJy and (i?i. 4 GHz) ^ 57%. 

If we imagine that our field had a sensitiv- 
ity (cti.4ghz) — 7.2 /iJy, six of our likely radio 
counterparts would fall below the 4.0 a limit of 
S1.4 GHz < 29yuJy and would not be detected. 
Four additional likely 20 cm counterparts would 
appear at the 4-5 a level and would be at high 
risk of not being detected due to the usual com- 
pleteness effects. Therefore, our radio counter- 
part identification rate would be 68+g%-78lg%. 
This range is only marginally greater than the av- 
erage value of 57%, and well within the scatter 
of the previous surveys. Therefore, we attribute 
our high radio identification rate to the extremely 
sensitive VLA map of this field, rather than to un- 
usual properties of 1.2 mm-selected sources at this 
depth. 

Because we expect ~ 5 spurious detections 
among our 41 sources with S/N > 4.0 and we 
find only 2-3 detections with unlikely or no ra- 
dio counterparts, there is little room left to ac- 
commodate a substantial, extremely high-redshift 



(z > 5) population of r adio-undetected SMGs (see 
also llvison et al.l 120051 ) . This work suggests that 
with a deep enough radio image, perhaps all SMGs 
might have their radio counterparts identified, au- 
guring well for upcoming deep surveys that exploit 
the dramatically expanded correlator bandwidth 
of the EVLA. 

We find that 7.3^11% (3/41) of our detec- 
tions have two likely radio counterparts (MAMBO 
ID # 3, 15, and 39). If we consider the fact that 
5% of all randomly chosen positions within our 
MAMBO map will have counterparts within 8" 
with P < 0.05, then we would expect to find 
a double radio counterpart rate of ~ 4.6% from 
chance associations. Previous studies have found 
that ~ 10% of SMGs host multiple likely radio 



eiy 

counterparts (se e, e.g., llvison et al.l 120021 12007 
Pope et al~ll2006f l. probably due to the effects of 



confusion within the submm/mm image, physi- 
cal interactions, or the extended jets of radio-loud 
AGN. Although our SMG sample in the LHN is 
too small to be able to constrain the fraction of 
multiple radio counterparts to better than ±5%, 
we note that the pair separations of the radio 
counterparts are 2.1", 7.7", and 7.4", and that 
two of the three MAMBO sources have deboosted 
flux densities in the top 25% of our sample. These 
results may be in agr e ement with the trend iden- 
tified in llvison et al.l (|2007l ) that multiple radio 
counterparts are preferentially associated with the 
brightest SMGs, and have pair angular separations 
AO ~ 2"-6". 

7.3. Redshift distribution 

As listed in Table [5] and detailed in Appendix 
IA1 of our 41 significant individual detections, two 
have o ptical spectroscopic redshi fts ( Polletta et al.l 
20061 : lOwen fc Morrisonl l2009ah. two have mid- 



IR spectroscopic redshifts (jFiolet et al. 2010h . and 
two have high-quality photometric redshifts based 
on Herschel far-IR photomet ry that we will denot e 



in what follows as "Zp hot " (jMagdis et al.l l2010h 



For those of the remaining 35 sources with robust 
or likely radio counterparts, we generally adopt 
the optical photometric redshift s (denoted z D hot 



i n wh at follows) determined by IStrazzullo et al 



( 2010 ) for the radio catalog of lOwen fc Morrison 



20081). The ex ception to this rule c omes for 



{■Zphot} to which Strazzullo et al. < 2010h assign a 
goodness-of-fit quality flag of "C" ; these redshifts 
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are less reliable, and in particular are more likely 
to manifest catastrophic errors. For such sources, 
as well as for the one 1.2 mm detection that lacks 
a radio counterpart altogether, we instead derive 
our own redshift estimates (z a ) using the radio- 
submillimeter spectr al index redshift indicator of 
Carilli fc Yunl (jl999l ): 



o 



350 



-0.24-[0.42 X (a ra dio - OWbmm) X log lf) (l + Z a ) 

, (16 ) 

(see also lCarilli fc Yunll2000l: lYun fe Carillill2002l) . 

For a r adio we use, in order of priority and availabil- 



ity, a 



90cm ^50cm 
20cm i u 20cmJ 



or 



-0.68. Weusealo 1 ™ 



only for 



sources with clean, unblended 50 cm detections. 
For these unblended 50 cm counterparts, we find 
an average value of (alo™) = —0.68 ± 0.06 (see 
Figure fT8|). in agreement with the average spec- 
tral index of SMGs in the LHE field of (a|g^) = 
-0.75 ± 0.06 (|lbar et all [20091 . l2010l) . We adopt 



this mean value of ct^ocm (—0.68) for redshift de- 
termination of sources with only 20 cm radio coun- 
terparts, or whose 50 cm counterparts are con- 
fused. For a S ubmm, we use the spectral index at 
1.2 mm of the fiducial high- redshift dusty galaxy 
SED (a su bmm = 3.2), as motivated in ^7.11 For 
our detection with no likely radio counterpart, we 
estimate a redshift lower bound by using the local 
£20 cm 4<7 upper limit in Equation 1161 

Figure [19] illustrates why we exclude C-quality 
photometric redshifts from our catalog. Plotted is 
the a\ 5 4 spectral index of the detections as a func- 
tion of Zphot- The points are co ded according to 



photo metric redshift quality flag (jStrazzullo et al 



2010). The points with the best photometric red 
shift fit quality (AA) are shown as black circles, 
followed by blue squares (A), green diamonds (B), 
and red triangles ( C). T he shaded region shows 
the Carilli fe Yunl ( 19991 ) relation for — a ra dio = 
0.52 — 0.80, where —0.52 represents the median 



a?,n ™ spectral index in the LHN field ([Owen et al 



2009 b]) and —0 .80 is the fiducial synchrotron value 
(|Condonlll992l ). The over-plotted lines show the 
empirical relations recovered by rcdshifting the 
SEDs of nearby star-forming galaxies M82 an d 



Arp 220 ({Klein et al.l 119881 : IScoville et all Il99lh 



The highest quality photometric redshifts agree 
with their galaxies' spectral indices in that they 
either follow the Carilli & Yun relation, or are con- 
sistent with an M82 or Arp 220 SED. In contrast, 
the C-quality photometric redshifts are scattered 
almost uniformly in z for a given af 5 4, demon- 



strating their lack of robustness. 

Figure [20] shows the redshift distribution of our 
catalog, including all spectroscopic, photometric, 
and a\ ^-estimated redshifts. It is apparent in Fig- 
ure [20] that the 17 spectroscopic and high-quality 
(Herschel and AA/A/B-grade optical) photomet- 
ric redshifts are biased towards lower redshifts. 
The median redshift for this 41% of our sample 
is ^median = 2.26, with an inter-quartile range of 
1.72-2.90. For all galaxies, z mc dian = 2.90, with 
an inter-quartile range of 2.33-3.70. This sys- 
tematic bias has two causes. First, the highest- 
redshift galaxies have the faintest counterparts, 
and will necessarily be detected in fewer optical 
bands, which results in a poorer fit. This trend 
is in contrast to the full radio catalog, for which 
the median z prio t is ~ 1 and the fraction of photo- 
metric redshifts with AA/A/B quality (~ 85%) 
is much higher than for our MAMBO s o urces . 
Second, the SEDs in the IStrazzullo et al.l (j2010l ) 
galaxy template library are most representative of 
nearby galaxies, potentially resulting in a poor fit 
if they are applied to high-z galaxies whose SEDs 
are not included in that library. 

The median redshift for our sample (^median = 
2.90) is larger than the m edian redshift deter- 
mined bv lPooe et all (|2006l) for a complete sample 
of 850 /im-selected SMGs with spectroscopic red- 
shifts (^median = 2.0). Although our redshift dis- 
tribution has greater uncertainties because it relies 
heavily on photometric redshifts and spectral in- 
dex redshi ft estimates, i t is in agreement with the 



results of Chapin et al. ( 20091 ). who have shown 
that, with high statistical significance, galaxies in 
a sample selected at 1 . 1 mm are detected at higher 
redshift (^median = 2.7) than those selected at 
850 yum. Our median redshift is also greater than 
that of the sample of 68 galaxies selected at 870 /xm 
from the LABOCA survey of the ECDF-S. Us- 
ing 17-band optical t hrough mid-IR photometry, 
Wardlow et al.l (|201 lh find z modian = 2.2. 



7.4. Spatial correlation with 20 cm sources 

The results from the w{9) and nearest neigh- 
bor analyses ( £)5.8[) suggest that our sources are 
clustered to some degree. The spike in w{9) on 
~ 4' scales is an intriguing result that is consis- 
tent with the visual impression of Figure [2] (i.e., 
the southeastern overdensity and the central void) 
and hints at the existence of large-scale structure 
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(LSS) in this field. To investigate whether the spa- 
tial distribution of our detections traces LSS that 
can also be seen at other wavelengths, we com- 
pare our source positions to the distribution of ra- 
dio sources within the LHN. In order to compare 
our MAMBO sources to radio sources at compa- 
rable redshifts, we include only radio sources with 
1.5 < Zphot < 4.5 (excluding all that only have a 
C-quality 2 p hot)- Additionally, we only consider 
radio sources with sizes greater than 1.0". These 
larger sources will be preferentially gas-rich merg- 
ers with extended star formation or radio-loud 
AGN, which we would naively expect to trace en- 
vironmental overden sities on the basis of studies at 
lower redshift (e.g.. iHill k Lillvlll99lb iBestl [2004 ; 
Best et all [2003: 1 Wake et al.ll2008l) . From a prac- 



samplc of high- z 20 cm galaxies. 

7.5. Resolving the 1.2 mm CIB 

At 1.2 mm, the cosmic infrared b ackground 
(CIB ) intensity is l v ~ 15- 24 Jydeg" 2 (jPuget et al 



tical standpoint, they can also be detected over 
the full area of the MAMBO map, allowing for 
a fair comparison; sources with 20 cm sizes < 1", 
in contrast, tend to be fainter, and therefore have 
systematically lower surface densities farther from 
the center of the VLA map. 

In Figure EU we plot our 1.2 mm source po- 
sitions over a 3'-resolution smoothed surface 
density map of the 30 7 radio sources in the 
lOwen k Morrison! (|2008t) VLA catalog that sat- 
isfy our selection cuts. We find that the dis- 
tributions of the two populations agree quite 
well: (1) the density map shows a deficiency 
of radio sources at the location of our central 
MAMBO void, (2) every radio source density 
peak ((S20cm) > 0.5arcmin -2 ) is associated 
with at least one MAMBO detection, and (3) 
~ 66% of our MAMBO detections are located 
in regions with higher than average ((£20 cm) — 
0.34arcmin ) radio source density, whose area 
comprises only 45% of the total. This strik- 
ing agreement seems to argue for a real physical 
correlation between our MAMBO detections and 
20 cm radio galaxies at sim ilar redshifts. Following 
Austermann et al. ( 2009h . we have also compared 



our SMG catalog to an identical number of homo- 
geneously distributed, randomly chosen positions. 
The experiment confirms the spatial correlation 
at a confidence level of 90%. However, when us- 
ing random positions derived from our simulated 
maps that take into account the spatially varying 
sensitivity, we find that our MAMBO detections, 
although significantly correlated with each other, 
are not significantly spatially correlated with this 



119961 : iFixsen et all Il998t >. By adding up the 
deboosted flux densities of our detections with 
S/N > 4.0, we recover ~ 0.58 Jydeg -2 of the 
CIB, or about ~ 3%. Figure [5] shows that our 
best fitting Schechter function estimate of the 
differential number counts is entirely consistent 
with the intensity of the CIB, while the power- 
law model is only marginal ly compatib l e with it. 
The analysis performed bv lScott et al. ( 2010l ) on 
ASTE/AzTEC data in the GOODS-S field finds 
that the best-fit power law model from their P(D) 
analysis can account for the CIB, although only if 
they integrate the counts past the cutoff used for 
that analysis. They also find that their Schechter 
function model is incompatible with the CIB, and 
can only recover ~ 30% of the 1.1 mm background 
when integrated down to S' l lmm = OmJy. These 
results may be due to their a priori choice of a 
faint-end power-law index 5' = 1.0. We fit for this 
parameter directly and find that a steeper value 
of 5' ~ 1.86 is optimal for our data and produces 
enough faint sources to account fully for the CIB 
light at 1.2 mm. 

8. Conclusions 

We have presented a 566 arcmin 2 map of 
the Lockman Hole North field with an aver- 
age optimally filtered point source sensitivity 
~ 0.75mJybeam _1 . By making use of previ- 
ously developed and original techniques to handle 
chopped bolometer array data, along with P(D)- 
based number counts and clustering analyses, we 
have assembled a comprehensive picture of the 
1.2 mm sky. Our results provide valuable new 
constraints for models of the evolution of dusty 
starburst galaxies through cosmic time. 

We detect 41 1.2 mm sources at S/N > 4.0 in 
our final map. Of these 41 detections, 38 have ro- 
bust or likely (P < 0.05) 20 cm radio counterparts, 
and 20 have robust or likely counterparts at 24 /im. 
Based on Monte-Carlo simulations, we expect ~ 5 
of these detections to be spurious, and only ~ 2 
20 cm counterparts with P < 0.05 to be chance 
associations. This result gives our MAMBO/LHN 
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map the highest single-field SMG radio counter- 
part identification rate ever observed (93±£%), 
which we have shown can be explained entirely by 
the extraordinary depth of our 20 cm VLA map. 
The enhanced sensitivity of the EVLA will be able 
to make high counterpart fractions routine for fu- 
ture SMG samples. Based on the spectroscopic, 
photometric, and radio/far-infrared spectral index 
redshifts of these counterparts, the median red- 
shift of our sample is z roe dian = 2.9, higher than 
has been determined for 850 /xm-selected SMG 
samples in fields with shallower VLA coverage. 

We estimate the number counts of 1.2 mm 
sources both directly and by using a P(D) anal- 
ysis and find a similar slope but a lower overall 
normalization relative to previous MAMBO sur- 
veys. However, our results are in close agreement, 
after a scaling in flux density, with those of recent 
surveys at 1.1mm. The compatibility of our di- 
rectly calculated counts and P(D) analysis with 
the constraint of the 1.2 mm CIB demonstrate the 
robustness of our results. In particular, we find 
that for 51.2 mm ^ 0.05 mJy the SMG differential 
number counts cannot keep rising with the faint- 
end slope observed for £1.2 nun > 0.05 mJy, and 
that the bright SMG population contributes at 
most a small fraction to the 1.2 mm CIB. 

The high resolution afforded by the IRAM 30 m 
telescope, the large extent of our map, and the use 
of analysis methods that thoroughly take into ac- 
count the negative residuals of the chopped triple- 
beam PSF have allowed us to demonstrate possi- 
ble clustering in the 1.2 mm population. The SMG 
correlation function, a nearest neighbors analysis, 
and, to a lesser extent, the spatial correlation of 
our significant detections with large radio sources 
over the same redshift range all suggest that our 
sample traces some degree of large scale structure 
at high redshift. Our work prepares the 1.2 mm 
waveband for the ALMA era by creating a bet- 
ter understanding of this population's statistical 
properties and setting new 1.2 mm constraints for 
galaxy evolution models. 
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A. Notes on individual detections 



A.l. MM J104700. 1+590109 = ID # 1 

Polletta et al. (2006) report an optical 



spec 



2.562 for this source. It also has a 70 /im counterpart with 



SVo nm = 10.4 ± 1.7 mJy at a distanc e of 3.2'' 
distance of 1.7" (jOwen et al.ll2011bh . 



and a 160 /um counterpart with S 



160 /im 



= 24.1 ± 1.9 mJy at a 



A. 2. MM J104631.4+585056 = ID # 3 



In addition to the robust 20 cm counterpart listed in Table [2 this source has an additional likely counter- 
part with 520cm = 30/zJy, separation 2.6", and P = 0.034, with which it is nearly blended. Neither radio 
counterpart has an estimated photometric redshift. Figure [2"2l shows 1.2 mm contours overlaid on a 20 cm 
cutout image that includes both counterparts. 

A.3. MM J104638.4+585613 = ID # 6 



(|2010h . from which it 



is 



We identify this source with LHN8 in the Herschel catalog of iMagdis et al 
se parated by 2 . 0" (P = 0.0084), and with SWIRE4.J10463 8.68+585612 . 5 = ID # L14 in the Spitzer sample 
of Fiolet et al. I fi2010h . who report a mid-IR z spcc = 2.03. iFiolet et al.l (|2009( ) report an on-off flux density 
measurement of Si. 2mm = 2.13 ± 0.71 mJy, which is consistent with our (non-deboosted) Si. 2mm = 2.7 ± 
0.5 mJy within the uncertainties. 

A.4. MM J104704. 9+585008 = ID # 9 



The radio counterpart to this source (P 
because it has a S/N ratio of 4.8. 



0.0043) is not in the catalog of I Owen & Morrisonl (|2008l ) 



A.5. MM J104556. 5+585317 = ID # 11 

We identify this source with LHN1 in the Herschel catalog of IMagdis et al. (2010). from which it is sepa- 
rated by 3.4 " ^ (P = 0.0057), and with SWIRE4.J104556.90+585318.8 = ID # Lll in the Spitzer sample of 
who rep o rt a m id -IR z Rpf , :C = 1 .95 in good agreement with the optical z p \ wt = 1-80 



P = 

Fiolet et all (l2010h. 



reported by IStrazzullo et al.l 

Si. 2mm 



(120101) . IFiolet et al.l (2009) report an on-off flux density measurement of 
3.08 ± 0.58 mJy, which is consistent with our ( non-deboosted) Si .2 mm = 3.4 ± 0.6 mJy within the 
uncertainties. This source also has a 160 /zm counterpart ( Owen et al.ll2011a ) with Sieo^m = 11.8 ± 1.5 mJy 
at a distance of 4.3" 



A. 6. MM J104728. 3+585213 = ID # 15 

This source has two radio counterparts; the primary counterpart listed in Table [2] has z p hot = 2.76, while 
the second (with P — 0.032) has z p hot = 1.06. In the 20 cm map, we see a quadruple radio source (see Figure 
I23| . At a low level of significance, the 1.2 mm emission ap pears to be elongate d in the same direction as the 
radio source(s). This system also has a 160 /im counterpart ( Owen et al.ll2011~bh with Sieo^m = 12.7+1.5 mJy 
at a distance of 6.8". 

A.7. MM J104610. 4+590242 = ID # 17 

This source is s epara ted by only ~ 15" from MM J104611. 9+590231 = ID # 39. The catalog of 
Owen fc Morrisonl (|2008l) includes a single radio source whose nominal position is midway between two 
very faint sources (one resolved, one unresolved) that are visible in the original 20 cm map (see Figure [24| . 
These two sources, with peak flux densities of 15.6 [iJy and 13.5 /zJy, are only identifiable because they lie 
very close to the center of the VLA map, where the local RMS noise is only 2.9 fiJy beam -1 . After attribut- 
ing the flux of the catalogued source to its two constituents, we find that one is a robust radio counterpart 
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for the MAMBO source (1.0" separation with P 
(3.3" separation with P = 0.077). 

A.8. MM J104617.0+585444 = ID # 20 



0.0099), while the other is probably a chance association 



This source has an unlikely 20 cm radio counterpart (P = .12) with S/N = 4.2 (5 2 o C m = 15.0 ± 3.6 £tJy), 
and is therefore not in the catalog of Owen fc Morrison! ( 2008 ). which includes only sources with S/N > 5.0. 



A.9. MM J104522. 8+585558 = ID # 26 

This source has no radio counterpart within 8" and no likely 24/im counterpart, although it does have an 
X-ray counterpart (see §6.5|) . 



A.10. MM J104620. 9+585434 = ID # 28 

This source has a likely 20cm ra dio counterpart (P = 0.02 9) with S/N = 4.3 (S 2 0cm = 24.8 ± 5.8 fiJy), 
and is therefore not in the catalog of Owen k, Morrison! ( 20081 ). which includes only sources with S/N > 5.0. 



A.ll. MM J104556. 1+590914 = ID # 29 



We i dentify this source with SDSS J104555. 49+590915. 9, an optically bright galaxy for which lOwen fc Morrison 



(|2009af ) report an optical z s 



= 0.044. Figure [25] shows a red optical image overlaid with 1.2 mm con- 
tours, which at a low level of significance are elongated in the same direction as the galaxy's stars. This 
~ 20 " source is heavily resolved at 20 cm and 50 cm. It is also detected at 160 fxra (lOwen et al. 2011b ) with 
Sieo^m = 15.8 ± 1.4mJy (at a separation of 6.6"), and at 250 /im (Herschel /SPIRE; ISmith et al.ll201l[) with 
£250 /mi = 133 ± 8mJy (at a separation of 3.6"). 

A.12. MM J104539.6+585419 — ID # 32 



We identify this source with LHN3 in the Hers chel catalog of Magdis et al. ( 2010h . from which it is 



separated by 5.2" (P = 0.016). iMagdis et all (j2010l ) estimate 



' — ' ' ppi^ 

SPIRE photometry, which we list in Table [2] rather than the IStrazzullo et al.l (J2.010) optical z p hot 

due to the close connection between far-IR and millimeter emission. 



2.40 on the basis of their PACS and 

1.32, 



A.13. MM J104608. 1+590744 = ID # 36 

This source has a 20cm ra dio counterpart (P = 0.04 1) with S/N = 4.9 (S20 cm = 16.1 ± 3.3 /iJy), and is 
therefore not in the catalog of lOwen fc Morrison (|2008h . which includes only sources with S/N > 5.0. 



A.14. MM J104610.8+585242 — ID # 37 

We identify this source with LHN4 in the Hers chel catalog of Magdis et al. ( 2010t ). from which it is 
separated by 2.4" (P = 0.013). IMagdis et all (|2010l) estimate z' phot = 1.72 on the basis of their PACS and 
SP IRE photometry , which we list in Table [2j this is in good agreement with the optical z p hot = 1-66 reported 
bv lStrazzullo et all <|2010h . 



A.15. MM J104611. 9+590231 — ID # 39 

This source is separated by only ~ 15" from MM J104610.4+590242 = ID # 17. It also has a pair of 
likely radio counterparts within an 8" search radius. Figur e I24 | shows 1.2 mm contou rs overlaid on a 20 cm 
cutout image. This source is not identified in the catalog of lOwen fc Morrisonl (|2008l ). 
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Fig. 1. — Weight map. Contours denote post- 
filtering RMS noise levels of 0.8, 1.0, 1.2, and 
1.4 mJy beam -1 in the "best" data (solid con- 
tours) and the "full" data (dotted contours). 
Thick dashed contour shows the map region used 
for the P(D) analysis of the "best" data. Cir- 
cles show the locations of detected sources with 
S/N > 4.0. The thick outer edge shows the ex- 
tent of the full map. The effective areas com- 
prising the "best" and "full" maps (RMS noise 
< 1.5 mJy beam -1 after filtering) are 514 and 
566arcmin 2 , respectively. 




Fig. 2. — The complete, optimally filtered S/N 
map of the "full" dataset, with white circles show- 
ing the locations of our 41 detections with S/N > 
4.0. 





Fig. 3. — Pixel flux distributions of the S/N maps. 
Left: Histograms are shown for both shuffled 
(thin) and jackknifed (thick) noise maps. Right: 
S/N map histogram for the "best" map (thick) and 
the "full" map (thin). Over-plotted in both panels 
is a Gaussian function with unit standard devia- 
tion and zero mean, normalized to the area under 
the histograms. All maps have a mean value con- 
sistent with zero (shown as the vertical line seg- 
ment), enforced by the chopped observing mode 
of the IRAM 30 m telescope and SAA reconstruc- 
tion. The "full" map has more pixels with high 
S/N because it reaches a higher sensitivity. The 
histograms were created with maps trimmed to a 
noise level of 1.5 mJy beam -1 . 
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Fig. 4. — The effective PSF (including nega- 
tive residuals) as a function of position across 
our "best" map. Red/blue represents a posi- 
tive/negative signal response. The thick solid con- 
tour encloses the area where the "best" map has 
RMS noise < 1.5 mJy beam -1 . 
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Fig. 5. — Pixel flux distribution showing the agree- 
ment between the data and the average best-fit 
(Schechter function) model from the P(D) anal- 
ysis. Points represent the PFD of the optimally 
filtered "best" data. The solid histogram shows 
the mean PFD of 100 random sky realizations of 
the best fitting number counts embedded in ran- 
dom jackknifed noise maps. The thin curve rep- 
resents the mean histogram of only the jackknifed 
noise maps. The single larger flux density bin that 
spanned 4-5 mJy for the P(D) analysis (see SJ3|) is 
shown here broken down into small bins matching 
the rest of the distribution. 
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Fig. 6. — Likelihood parameter spaces for dif- 
ferent parametrizations of the 1.2 mm differential 
number counts. The points mark the best fit- 
ting parameters and the marginalized 68% double- 
sided error bars. The contours bound the 68% 
and 95% confidence regions found through Monte 
Carlo simulations. The shaded bands show the 
regions in parameter space that reproduce the ob- 
served intensity of the 1.2 mm CIB. assuming a 
lower-limit flux density cutoff of 0.05 mJy. Left: 
Power law model, with best fitting parameters 
A^mjy = 19-7±|;| dcg^mJy- 1 and S = 3.14^. 
Right: Schechter function model, with best fit- 
ting parameters N± mJy = 14. 5^7 deg _2 mJy _1 
and 5' = 1.861^2. 
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Fig. 7. — Average angular separation between in- 
jected and recovered source locations as a func- 
tion of flux density. The injection process was 
performed 1000 times for each flux density value. 
The scatter increases at low flux densities because 
of incompleteness. 
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Arc Minutes 

Center: R.A. 10 46 07.43 Dec +59 01 16.6 

Fig. 8. — The change in S/N of our 17 compar- 
ison sources in the "best" and "full" datasets. 
Filled circles show eight sources in the "best" 
map with S/N > 4.5. In the "full" map, 
we recover all of these detections, plus oth- 
ers shown as empty circles, above a threshold 
S/N > 5.0. The solid/dashed contour denotes the 
lmJybeam -1 RMS noise threshold in the maps 
of the "best" / "full" data. The numbers near each 
detection show the S/N ratio of that source in the 
"best" data — > "full" data. 
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Fig. 9. — Flux density comparison between the 17 
highest S/N detections in the "full" data and their 
counterparts in the "best" data (these sources are 
plotted in Figure [5]). The solid lines show the best 
fit slope with 2cr uncertainties of a line constrained 
to cross the origin. For illustration, the dotted line 
has unit slope. The chi-square minimization with 
x and y erro rs was performed u sing the IDL script 



20081 ). The best fitting 



mpfit.pro (jMarkwardt 
slope is Tji = 1.06 ± 0.07 (±lcr uncertainties). 
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2 

S, [mjy] 

(On-Off Photometry, Fiolet et 01. 2009) 

Fig. 10. — Flux de nsity comparison b etween on-off 
photometry- mode ( Fiolet et al. 20091 ) and on-the- 
fly mapping (this work) of 13 Spitzer-selected star- 
burst galaxies in the LHN that lie within our "full" 
map coverage. Symbols with thick (thin) lines 
repres ent 2 (11) of the sources from iFiolet et all 



(l2009h with significant detections (non-detections) 
in our full map. 
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Fig. 11. — Integral numbers of spurious detections 
function of S/N within jackknifed noise maps, 
shuffled noise maps, and maps of random Gaussian 
numbers. We expect to detect 0.8 spurious sources 
with S/N > 4.5 and ~ 5.4 spurious sources with 
S/N > 4.0. The curve shows the expected num- 
ber of excursions above a given S/N level within 
a Gaussian r andom field in the limit of high S/N 
(|Adler Ill98lh . 
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Fig. 12. — Completeness as a function of flux den- 
sity. The circles/diamonds represent 4.5cr/4.0er 
source extraction thresholds for detection. The 
number count calculation uses a threshold of 4.5a, 
while our source catalog includes sources down 
to 4.0cr. Horizontal lines represent 80% and 50% 
completeness limits. 
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Fig. 13. — Posterior probability distributions of 
four characteristic S/N regimes. The dotted lines 
represent the Gaussian probability distributions 
for each of the four measurements (the likeli- 
hoods), all assuming a — 0.6 mjy and varying 
mean. The dashed lines represent the prior flux 
distribution constructed from Monte Carlo simu- 
lations, and the solid lines represent the the nor- 
malized posterior probability distributions. The 
brackets above the curves denote the peak likeli- 
hood values (the deboosted flux densities) along 
with the left and right 38% confidence intervals. 
When the S/N is too low for either the left or the 
right confidence interval to converge, as is the case 
in the final panel, we use instead the analytic for- 
mula from Equation 1121 
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Flux Density [mjy] 



Fig. 14. — Monte Carlo simulation to test our 
Bayesian flux boosting correction. The histogram 
shows the average number counts model used to 
inject sources into realistic noise maps. The circles 
show the average number counts calculated from 
the raw counts by using the Bayesian method to 
deboost the flux density of each recovered source. 
Error bars repr esent only the Poisson uncertainty 
in the average (lGehrelslll986l) . 



Fig. 15. — Left: 1.2 mm integral number counts, 
compared to the P(D) analysis best fit models 
and other observed number counts from the lit- 
erature. Filled black circles show our Bayesian 
deboosted number counts with 95% confidence er- 
ror bars including the Poisson uncertainty and 
the uncertainty in the completeness correction. 
The solid filled region shows the 95% confidence 
region for the best fitting Schechter function 
from the P(D) analysis. The cross-hatched re- 
gion shows the same for the best fitting power 
law. MAMBO counts: LH and ELAIS-N2 
(gree n squares and blue triangles 



Greveetal 



2004h . COSMOS (green circles; iBertoldi eta] 
20071 ). LA BOCA counts: ECDF-S (purple stars; 



Weifi et al.l l 2009h. AzTE C counts: GOODS-S (or- 
ange ci rcles: IScott et all |2010|). SHADES (yellow 
circles; lAustermann et al.l l2010l). COSMOS (red 
circles; Iseott et al.ll2008l ). SCU BA le nsing clus- 



ter counts from Knudsen et al. (120081) are shown 
as red stars. The counts at 850 yum (SCUBA), 
870 /mi (LABOCA), and 1.1mm (AzTEC) have 
been scaled to 1.2mm (see §7.ip . Right: 1.2mm 
differential number counts and P(D) models com- 
pared to the predictions of galaxy evolution mod- 
els. Lines re present different d i fferen tial counts 
predicti ons by Rowan- Robinson! ( 2009f ). blue dot- 
dashed; Valiante et al. ( 20091 ), gree n short-dashed 



Bet hermin et al.l (|201l[ ). red solid; iMarsden et al 



(|201ll ) , cyan long-dashed. Models with predictions 
only for the 1 . 1 mm waveband were scaled in order 
to compare to our observations (see £|T. ID . 
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Fig. 16. — Landy-Szalay correlation function es- 
timator w(6) as a function of angular separation. 
Top panel uses our 41 detections with S/N > 4.0. 
Bottom panel uses 27 sources with S/N > 4.5. 
The vertical dashed line shows upper limit in 9 
on the clustering suppression effect introduced by 
chopping (see ^5.8|) . Open squares show the re- 
sults from using random positions to check the 
zero-clustering baseline; open circles show the raw 
clustering signal, which are corrected for the zero- 
clustering baseline to deliver the filled circles. 
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Fig. 17. — Integral distribution of nearest neigh- 
bors. Solid line shows the distribution from 
Monte-Carlo generated random positions. Circles 
show the distribution of our 41 significant detec- 
tions. The vertical line segment denotes the max- 
imum difference between the two distributions. A 
Kolmogorov-Smirnov test rules out the hypothesis 
that our sources are drawn from a random distri- 
bution at 95% confidence. 
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Fig. 18. — 50/20 cm spectral index vs. 20 cm flux 
density for SMG radio counterparts. Solid circles 
represent the 20 unblended 50 cm radio counter- 
parts to our MAMBO detections. Open circles 
(squares) represent non AGN-domina ted (AGN- 



dominated) SMGs in the LHE field (jlbar et al 
l2010h . 



33 



15 



io 10 ■ 



ir 0.6 




Fig. 19.— 350/1.4 GHz spectral index (a? 5 4 °) vs. 
z. Black circles, blue squares, green diamonds, 
and red triangles have z p hot with ht qualities 
of AA, A, B, and C, respectively. The shaded 
area shows th e expe cted behavior of a\ 5 ® for the 
Carilli fc Yunl (|l999h spectral index redshift indi- 
cator, with — a ra dio ranging from 0.52-0.80. The 
solid and dashed lines show the empirical relations 
obtained by redshifting the SEDs of nearby star- 
burst galaxies Arp 220 and M82, respectively. 
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Fig. 20. — Redshift distribution of significant de- 
tections. The light grey histogram includes all 
sources. The dark grey histogram includes only 
sources with spectroscopic or high-quality (Her- 
schel or AA/A/B-grade optical) photometric red- 
shifts. 




Arc Minutes 

Center: R.A. 10 46 07.43 Dec +59 01 16.6 

Fig. 21. — Distributions of individually detected 
MAMBO sources and the smoothed surface den- 
sity of 20 cm radio sources with 1.5 < z p hot < 4.5 
(excluding radio sources with C-quality photomet- 
ric redshifts) and sizes > 1". The area shown is 
the region in the "full" map used for source ex- 
traction. 
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a (J2000) 

Fig. 22. — 1.2 mm contours (multiples 

of 1.0 mJy beam -1 ) for MAMBO source 
MMJ104631.4+585056 = ID # 03, overlaid 
on 20 cm map showing double counterpart. The 
black circle and small filled ellipse at lower 
right represent the 15.6" MAMBO beam (after 
filtering) and the 1.63" x 1.57" VLA beam, 
respectively. 




a (J2000) 

Fig. 23. — 1.2 mm contours (multiples 

of 1.0 mJy beam -1 ) for MAMBO source 
MMJ104728.3+585213 = ID # 15, overlaid 
on 20 cm map showing multiple counterparts. 
Other notation is as in Figure [2U 




a (J2000) 

Fig. 24— Nearly blended MAMBO sources 
MMJ104610.4+590242 (ID # 17), which 
has one likely 20 cm counterpart, and 
MMJ104611.9+590231 (ID # 39), which has 
two likely 20 cm counterparts. 1.2 mm contours 
(multiples of 0.7 mJy beam -1 ) are overlaid on 
20 cm greyscale. White circles are centered on 
the positions of the extracted MAMBO sources; 
other notation is as in Figure H2J 



35 




59 08 45 - 



10°46 m 00 



58 s 56 s 

a (J2000) 



Fig. 25. — 1.2 mm contours (multiples 

of 0.8 mJy beam" 1 ) for MAMBO source 
MM J104556. 1+590914 = ID # 29, overlaid on a 
DSS red image of SPSS J104555 49+590 915.9 at 
z spcc = 0.044 (jOwen fc Morrisonl E)09al ) . Other 
notation is as in Figure [22] 
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Table 1 
mambo observations 









Tracking 


Bad refraction 


Chop throws 


Season 


Maps 


Hours 


jitter? 


correction? 


(arcsec) 


2006 winter 


76 


78.9 


X 


X 


36/48 


2008 summer 


12 


12.3 


X 




42/36 


2008 winter 


39 


40.3 






42/36 


2009 summer 


8 


8.0 






42/36 


2009 winter 


52 


53.0 






42/36 


Total "best" 


99 


101.3 








Total "full" 


183 


192.5 









Note. — "Best" data use maps from winter 2008, summer 2009, and winter 
2010. "Full" data use maps from all seasons in the table. 
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Table 2 
mambo detections 











^iDcboosted 






^20 




f 


°1Q cm 


h 

•^50 cm 


•^90 cm 




k 

•^24 ^m 


ID 


Source name 


(mJy) 


(mJy) 


(mjy) 


p« o) d 


z 


(arcsec) 


-^20 cm 


C/^Jy) 


(MJy) 


(^Jy) 


j 

-^24 ^m 


(MJy) 


S/N > 4.5 


1 


MM J104700. 1+590109 


3.7+0.8 


4.1+0.6 


o r + O-ti 
°°-0.6 


< 0.01 


2.562 1 


2 





0.0025 


278+6 


490 ± 11 


687+72 


0.0018 


1280 ± 9 


2 


MM J104627. 1+590546 


4.5+0.8 


4.7+0.7 


, q + 0.7 
6 -°-0.T 


< 0.01 


4.29 m 


3 





0.037 


35+6 


49 ± 23* 




0.11 


41 + 6 


3 


MM J104631. 4+585056 


6.1+1.8 


4.7+0.7 


Q q + 0.8 

s -°-o.r 


< 0.01 


1.8" 





.5 


0.0001' 


458+9 


937 ± 10 


1633+85 






4 


MM J104607. 4+585413 


2.8+0.8 


3.2+0.5 


9 7 +0.5 
Z -' -0.5 


< 0.01 


4.4" 





.5 


0.0025 


29+3 


50 ± 11 








5 


MM J104725. 2+590339 


4.9+0.9 


5.2+0.8 


4 +0 ' 8 


< 0.01 


3.00'" 


3 


4 


0.028 


51 + 5 


115 ± 10 




0.0099 


395 ± 26 


6 


MM J104638. 4+585613 


3.1+0.7 


2.7+0.5 




< 0.01 


2.03° 


2 


1 


0.0043 


159+5 


321 ± 9 


442+70 


0.0074 


662 ± 8 


7 


MM J104700. 1+585439 


2.8+0.7 


2.8+0.5 


9 , + 0.4 
z -°-0.5 


< 0.01 


4.2" 


2 





0.016 


41+4 


60 ± 10 




0.10 


329 ± 18 


8 


MM J104633. 1+585159 


4.5+1.2 


3.4+0.6 


9 7 +0.6 
z - ' -0.6 


< 0.01 


3.3" 





4 


0.0003 


97+9 


145 ± 10 




0.0022 


342 ± 8 


9 


MM J104704. 9+585008 


5.6+1.5 


5.1+0.9 




< 0.01 


3.9" 





6 


0.0043 s 


23+5 


67+ 11 








10 


MM J104622. 9+585933 


3.6+0.7 


2.9+0.5 


2 4+ ' 5 

Z - -0.5 


< 0.01 


2.6" 


2 


.1 


0.0083 


78+5 


153 ± 11 


383+73 


0.14 


221 ± 8 


11 


MM J104556. 5+585317 


3.5+0.9 


3.4+0.6 


n 7+O.6 
z - 1 -0.6 


< 0.01 


1.95° 


3 


.5 


0.0057 


314+10 


427 ± 9 


662+74 


0.017 


684 ± 7 


12 


MM J104448. 0+590036 


5.1+0.9 


3.5+0.6 


9 7+O.6 
z - 1 -0.7 


< 0.01 


2.16 m 


3 





0.0049 


273+13 


421 ± 20 


815 + 140 


0.037 


597 ± 27 


13 


MM J104609. 0+585826 


2.7+0.7 


2.7+0.5 


9 1 +0.5 
^■ 1 -0.5 


< 0.01 


1.14'" 


3 


.1 


0.0067 


197+3 


461 ± 10 


1101 + 151 


0.28 


128 ± 8 


14 


MM J104636. 1+590749 


4.3+0.8 


4.3+0.8 




< 0.01 


2.26*" 


2 


4 


0.0079 


97+3 


215 ± 14 




0.0055 


596 ± 7 


15 


MM J104728. 3+585213 


5.1 + 1.1 


4.5+0.9 




< 0.01 


2.76 m 


2 


4 


0.0044' 


180+15 


245 ± 12 




0.0081 


834 ± 23 


16 


MM J104720. 9+585151 


2.3+1.0 


3.8+0.7 


9 7 +o.s 

z - ' -0.8 


< 0.01 


4.9" 


1. 


5 


0.0088 


47+6 


55 ± 15 




0.13 


307 ± 9 


17 


MM J104610. 4+590242 


3.4+0.7 


2.6+0.5 


a.o±S;» 


< 0.01 


4.0" 


1. 


1 


0.0099 s 


27+3 


< 45 




0.025 


293 ± 19 


18 


MM J104655. 7+585000 




4.6+0.9 




0.017 


1.30'" 


5. 





0.024 


104+6 


202 ± 11* 


385+73 


0.026 


679 ± 9 


19 


MM J104502. 1+590404 




2.6+0.5 


2 +0 (3 


< 0.01 


4.1" 


1. 


.8 


0.020 


28+4 


47 ± 12 








20 


MM J104617.0+585444 




2.3+0.5 


i 7+0.5 
1 - 1 -0.5 


< 0.01 


> 4.6" 


5. 





0.12 s 


15 + 4 


< 33 




0.027 


67 + 8 


21 


MM J104530. 3+590636 




2.4+0.5 


1 Q + 0.5 


< 0.01 


3.1" 








0.0057 


36+3 


82 ± 12 




0.017 


196 ± 8 


22 


MM J104603. 8+590448 




2.7+0.6 


2 +0 (3 


< 0.01 


1.44'" 


3 


4 


0.0094 


165+5 


355 ± 15 


485+72 


0.019 


595 ± 9 


2:5 


MM J104641. 0+585324 




2.2+0.5 


i 7+0.5 
i ' ' -0.5 


< 0.01 


3.6" 


1. 


6 


0.016 


31+4 


56 ± 12 








24 


MM J104500. 5+590731 




2.4+0.5 


1 s + 0.6 


< 0.01 


3.24'" 


2 


3 


0.011 


67+7 


132 ± 10 




0.089 


264 ± 18 


25 


MM J104540. 3+590347 




2.3+0.5 


i 7+O.6 
i ' ' -0.6 


< 0.01 


3.5" 


2 


3 


0.036 


24+6 


53 ± 15 




0.0071 


119 ± 7 


26 


MM J104522. 8+585558 




2.5+0.6 


1 s + - 7 


< 0.01 


> 5.0" 








< 14 


< 40 




0.11 


239 ± 20 


27 


MM J104702. 4+585102 




3.1+0.7 


2 +tKS 


0.018 


2.9" 





2 


0.0002 


77+12 


49 ± 12 




0.030 


237 + 9 


4.0 < S/N < 4.5 


28 


MM J104620. 9+585434 




2.1+0.5 


i r + 0.5 
J -°-0.5 


< 0.01 


3.8" 


2 


1 


0.029 s 


25 + 6 


< 60 








20 


MM J104556. 1+590914 




2.8+0.7 


i s+o.s 


0.022 


0.044 p 


4 


.7 


0.0092 


307+39 


567 ± 85 




0.0012 


4838 ± 23 


30 


MM J104510. 3+590408 




2.2+0.5 


1 c + 0.6 
1 -°-0.6 


0.011 


2.4" 


1. 


2 


0.0034 


74+6 


155 ± 10 




0.012 


56 + 7 


31 


MM ,1104624.7+585344 




2.1+0.5 


1 c + 0.6 
J -°-0.5 


< 0.01 


2.90 m 


2 





0.027 


43+3 


72 + 9 




0.079 


236 ± 7 


32 


MM J104539. 6+585419 




2.5+0.6 


i 7+O.8 
1 - 1 -0.7 


0.023 


2.40" 


5. 


2 


0.056 


46+3 


113 ± 9 




0.059 


485 ± 7 


33 


MM J104535. 5+585044 




5.2+1.2 


2.71- 


0.20 


3.7" 





3 


0.049 


70+9 


101 ± 10 




0.19 


245 ± 7 


34 


MM J104453. 7+585838 




2.4+0.6 


i-5i° ; 7 7 


0.022 


3.5" 


3 





0.031 


40+8 


65 ± 10 




0.13 


76 + 7 


35 


MM J104717.9+585523 




2.4+0.6 


1 5+ - 7 
J -°-0.7 


0.025 


3.71 m 


2 


4 


0.014 


60+5 


68 ± 11 




0.099 


297 ± 23 


36 


MM J 104608. 1+590744 




2.9+0.7 




0.052 


4.5" 


2 


1 


0.041 s 


16 + 3 


< 47 




0.11 


396 ± 39 


37 


MM J104610. 8+585242 




2.3+0.6 


"IS:? 


0.030 


1.72" 


2 


3 


0.0052 


160+7 


320 ± 10 


450+75 


0.026 


401 ± 8 


38 


MM J104444. 5+590817 




2.8+0.7 


1 = + 0.9 
1 -°-0.9 


0.052 


3.6" 


2 


.8 


0.017 


68+16 


77 ± 10 




0.12 


153 ± 8 


39 


MM J104611. 9+590231 




2.1+0.5 




0.023 


2.6" 


2 


2 


0.024' 3 


34+3 


99 ± 13 




0.024 


634 ± 8 


40 


MM ,1104658.7+590633 




3.0+0.8 




0.074 


2.9" 


2 


6 


0.017 


56+4 


102 ± 11 




0.032 


229 ± 8 


41 


MM J104600. 7+585502 




2.1+0.5 




0.025 


2.6" 


3 


.1 


0.018 


71 + 11 


67 + 9 









Note. — P values in boldface type denote a likely counterpart (P < 0.05). 
a Raw flux density extracted from our "best" map (see EfB~2l . 
b Raw flux density extracted from our "full" map(scc *|5.2| >. 

c Flux density extracted from our "full" map, and corrected for flux boosting (sec 35. 6( . 
d Total probability that the deboosted flux density is < m Jy (sec jj5.6J 
c Angular separation between MAMBO detection and 20 cm counterpart. 
t Probability of chance association with 20 cm counterpart. 
s Flux density of 20cm counterpart (Owen Sz Morrison 2008). 

Flux density of 50 cm counterpart J Owen et al.|[20TTa|) . 
*Flux density of 90cm counterpart (Owen ct al. 2009b). 
J Probability of chance association with 24 fim counterpart. 
k Flux density of 24 fim counterpart i jOwen et al.H2011bl) . 
'zspec from lPollotta et alj ( 1200611 . 
m 2phot from Strazzullo ct al. (201(1). 

n z a estimated from[Carilli & Yun (199^) spectral index (sec (|7..1t . 
°z Bp ,c from lFiolet et al.l (120101) . 
p z spoc from lOwen fe Morrison! d2009aT l. 
q z phot from lMagdis et al.l d2010T) . 

r Thcrc are two radio counterparts within 8" of MAMBO source with P < 0.05; see Appendix [A] 
s The catalog of[Owcn fo Morrison! (1200811 does not contain this radio counterpart (sec 36. H . 
*Flux density is uncertain due to blending (sec £|6.2t . 



Table 3 
1.2 mm number counts 



Differential Integral 



Flux Bin 
(mJy) 


Flux Density 
(mJy) 


dN/dS 
(deg -2 mJy -1 ) 


Flux Density 
(mJy) 


N{> S) 
(deg- 2 ) 


1.68-2.14 
2.14-2.59 
2.59-3.05 
3.05-3.51 
3.51-3.97 


1.91 
2.36 
2.82 
3.28 
3.74 


435±11 
12811™ 

1561SJ5 1 
18+8° 

oz -45 


1.68 
2.14 
2.59 
3.05 
3.51 


3661 2 * 2 
166l^ 8 
1081°° 

07+58 

— 22 

28±ff 



Table 4 
Area surveyed 



Area Depth HPBW 

Reference Instrument Field (deg 2 ) (mjybeam^ 1 ) (arcsec) 



Greve et al. (20041 MAMBO 1.2mm ELAIS-N2 & LH 0.099 0.8 11 

Bertoldi et al. f20071 MAMBO 1.2 mm COSMOS 0.11 1.0 11 

Scott et al"72*008) AzTEC/JCMT 1.1 mm COSMOS 0.15 1.3 18 

Perera et al. (2008) AzTEC/JCMT 1.1 mm GOODS-N 0.068 0.96-1.16 18 

Weifi et al. (2009) LABOCA/APEX 870 fj.ni ECDF-S 0.25 1.2 19.2 

Austermann et al (2010) AzTEC/JCMT 1.1 mm LHE & SXDF 0.5 1 18 

Scott et al. (2010) AzTEC/ASTE 1.1 mm GOODS-S 0.14 0.48-0.73 30 

Hatsukadc el al. (2011) AzTEC/ASTE 1.1 mm AKARI, SXDF, fc SSA 0.25 0.32-0.71 30_ 

This Work MAMBO 1.2 mm LHN 0.16 0.75 11 



Table 5 

20 cm counterpart identification rates 









01.4GHz 




Reference 


Instrument 


Field 


(/iJy beam -1 ) 


ID rate 


Borvs et al. (2004) 


SCUBA 850 Aim 


GOODS-N 


9.0 


58% (11/19) 


Ivison et al. (2007) 


SCUBA 850 /xm 


LHE & SXDF 


4.2 & 7 


52% (62/120) 


Bertoldi et al. f2007) 


MAMBO 1.2 mm 


COSMOS 


8.5 


73% (11/15) 


Scott et al. (2008) 


AzTEC 1.1mm 


COSMOS 


10.5 


44% (12/47) 


ChaDin et al. C2009) 


AzTEC 1.1mm 


GOODS-N 


4.5 


76% (22/29) 


Bigas et al. f2011) 


LABOCA 870 /im 


ECDF-S 


6.5 


37% (47/126) 


This Work 


MAMBO 1.2 mm 


LHN 


2.7 


93% (38/41) 
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